Historic,  Archive  Document 

Do  not  assume  content  reflects  current 
scientific  knowledge,  policies,  or  practices. 


United  States 
Department  of 
Agriculture 

Forest  Service 

Intermountain 
Research  Station 

General  Technical 
Report  INT-208 

July  1986 


Ms 


Supplement  to  the 
User's  Guide  for  the 
Stand  Prognosis 
IVIodel— Version  5.0 


William  R.  Wykoff 


CR=fe[CCF"] 

B=^j^''xe-xy2dx 

ln[AHj=-f,[AD,H,D] 

■    MORT=|+exp  [-BiXi] 

fin  [bAl]  -fe  [D.b.h.,HabitaLt,Crown] 


C3 


CO 


THE  AUTHOR  CONTENTS 

WILLIAM  R.  WYKOFF  is  a  research  forester  with  the  Page 

Forestry  Sciences  Laboratory  in  Moscow,  ID.  He  has  Introduction   1 

worked  on  growth  and  yield  of  Northern  Rocky  Changes  in  Predictions  of  Tree  Development  ...  i 

Mountain  conifers  since  1974.  Much  of  this  work  has  Snnall-Tree  Models   2 

focused  on  development  and  implementation  of  the  Increment  for  Large  Trees   5 

Prognosis  Model.  Individual  Tree  Mortality  Models   9 

Random  Effects   11 

RESEARCH  SUMMARY  interpretation  of  Input  Stockability  Data  11 

Calibration  of  Growth  Models   11 

Version  4.0  of  the  Prognosis  Model  was  released  in  Attenuation  of  Model  Calibration   12 

September  1981.  Since  then,  a  regeneration  establish-  Assigning  Initial  Values  for  Missing  Heights  .  12 

ment  model  has  been  completed  and  small-tree  incre-  Computing  Scribner  Board  Foot  Volume  with 

ment  models  have  been  greatly  refined.  The  COVER  Variable  Merchantability  Standards   12 

model  has  also  been  added  to  predict  shrub  develop-  jree  Record  Compression  14 

ment  and  total  canopy  cover.  Thus,  the  representation  Changes  in  Features  that  Were  Part  of 

of  the  vegetative  component  of  the  stand  is  basically  Version  4.0  "  14 

complete  and  the  Stand  Prognosis  Model  can  be  Seeding  the  Random  Number  Generator  14 

linked  more  readily  to  models  for  nontimber  resources.  Generating  a  List  of  Projected  Tree  Records  .  14 

New  management  options  have  been  added  to  the  Additional  Flexibility  in  Field  Sampling 

system,  and  an  Event  Monitor  increases  the  flexibility  Procedures  15 

for  scheduling  management  activities.  A  compression  jree  Records  Can  Be  Completely  Removed 

or  classification  algorithm  enhances  program  effi-  in  Thinning  15 

ciency  by  combining  tree  records  that  are  similar  with  o.b.h.  Limits  Added  to  Thinning  Options  16 

regard  to  attributes  that  influence  growth  predictions.  Access  to  the  COVER  Extension  16 

Finally,  there  have  been  numerous  improvements  in  Interpreting  Stand  Information  Data  17 

the  biological  models.  Mortality  Multipliers  17 

This  report  is  a  supplement  to  the  user's  guide  for  xree  Records  17 

the  Stand  Prognosis  Model  (Wykoff  and  others  1982).  New  Options  18 

Options  that  were  available  in  version  4.0  may  still  be  Event  Monitor  18 

invoked  in  the  manner  described  in  the  user's  guide.  Variable  Merchantability  Standards  for  the 

New  options,  new  models,  and  modifications  to  exist-  Board  Foot  Volume  Equations  18 

ing  models  are  herein  described  as  incorporated  in  Descriptive  Statistics  for  Input  Data  19 

version  5.0.  Starting  Projections  from  Bare  Ground   19 

Compressing  the  Tree  List  19 

Simulating  Top  Mortality  20 

Changes  in  Program  Output  20 

Conclusions  21 

References  26 

Appendix  A:  Summary  of  Codes  Used  in  the 

Prognosis  Model  27 

Appendix  B:  Summary  of  Keywords,  Associated 

Parameters,  and  Default  Conditions  29 

Rules  for  Coding  Keyword  Records   29 

Controlling  Program  Execution   29 

Entering  Stand  and  Tree  Characteristics  30 

Modifying  Tree  Volume  Calculations   31 

Specifying  Management  Activities   32 

Controlling  Program  Output  33 

Linkage  to  Prognosis  Model  Extensions  34 

Growth  Prediction  Modifiers  and  Special 
Input/Output  Options  34 


July  1986 
Intermountain  Research  Station 
324  25th  Street 
Ogden,  UT  84401 


Supplement  to  the  User's 
Guide  for  the  Stand  Prognosis 
Model— Version  5.0 

William  R.  Wykoff 


INTRODUCTION 

Version  5.0  of  the  Prognosis  Model  for  Stand  Development  (Stage  1973)  was 
released  in  July  1984.  This  version  contains  many  enhancements  that  follow  the 
release  of  version  4.0  in  September  1981  (Wykoff  and  others  1982).  Biological 
models  were  reformulated  to  correct  apparent  deficiencies,  and  submodel 
parameters  were  reestimated  to  take  advantage  of  new  data.  Options  have  been 
added  to  increase  the  scope  of  the  program  and  ease  user  interaction.  In  addi- 
tion, several  new  extensions  have  been  developed. 

In  order  to  maintain  a  semblance  of  stability  for  program  users,  very  few  of 
these  changes  were  incorporated  into  version  4.0.  (Copies  of  version  4.0  source 
code  were  distributed  to  29  public  and  private  organizations.)  Programming  and 
model  formulation  errors  were  corrected  and  distributed,  but  version  4.0 
produces  essentially  the  same  results  today  that  it  did  when  first  released. 

Developments  from  the  period  1981  through  1984  have  been  implemented  in 
Prognosis  Model  version  5.0.  These  developments  include: 

-  Revision  of  the  small-tree  height  increment  model  and  the  large-tree  d.b.h. 
increment  model. 

-  Replacement  of  the  small-tree  d.b.h.  increment  model,  the  mortality  model, 
and  the  crown  dubbing  procedure. 

-  Replacement  of  the  procedure  for  calibrating  the  large-tree  d.b.h.  increment 
and  smzdl-tree  height  increment  models  so  that  the  average  predicted  increment 
matches  the  average  observed  increment,  regardless  of  sampling  design. 

-  Modifications  to  the  DESIGN,  and  THINning  options  that  increase  pro- 
gram flexibility. 

-  Introduction  of  new  options  such  as  NOTREES,  HTGSTOP,  TOPKILL, 
COMPRESS,  and  an  Event  Monitor. 

-  Linkages  to  a  western  spruce  budworm  model,  a  regeneration  establishment 
model,  COVER,  and  a  pau-allel  processor. 

This  supplement  to  the  User's  Guide  describes  differences  in  the  behavior  of 
biological  models,  provides  reference  information  on  using  new  extensions,  and 
describes  new  options  and  modifications  to  existing  options.  Benchmarking 
procedures  and  results  of  test  runs,  including  comparisons  with  version  4.0,  will 
be  reported  separately.  These  benchmarks  will  provide  a  mechanism  for  the  sys- 
tematic evaluation  of  future  model  changes. 

CHANGES  IN  PREDICTIONS  OF  TREE  DEVELOPMENT 

Version  5.0  retains  two  of  the  five  submodels  that  were  part  of  version  4.0 
(the  large-tree  height  increment  and  crown  change  models).  Three  submodels 
have  been  modified  or  replaced  (the  mortality  model,  the  small-tree  height  incre- 
ment model,  and  the  large-tree  d.b.h.  increment  model).  In  addition,  new  models 
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for  predicting  d.b.h.  increment  and  crown  change  in  small  trees  (d.b.h.  <  3.0 
inches)  were  implemented.  New  models  and  modifications  are  briefly  described 
below.  Coefficients  are  listed  for  the  models  that  differ  from  version  4.0 
formulations. 

Small-Tree  Models  Small-tree  growth  has  been  modified  considerably  by  changes  in  the  mechan- 

isms for  prediction  of  height  increment,  d.b.h.  increment,  and  change  in  crown 
ratio.  The  small-tree  height  increment  model  (eq.  1)  now  includes  relative  tree 
size  effects  and  has  been  expanded  to  apply  to  additional  habitat  types.  These 
effects  were  borrowed  from  the  large-tree  d.b.h.  increment  model  that  is  part  of 
version  4.0  (Wykoff  and  others  1982).  For  trees  with  d.b.h.  less  than  2  inches, 
the  height  increment  prediction  is  based  entirely  on  equation  1.  For  trees 
between  2  and  10  inches  d.b.h.,  small-  and  large-tree  height  increment  predic- 
tions are  averaged  as  described  for  version  4.0. 

In  (HTG2)  =  HAB  +  LOC  +  0.22157-SL.cos(ASP)  -  0.12432-SL.sin(ASP) 

-  0.10987-SL  +  bi-ha(HT)  +  bg-CCF  +  bg-IBAL/lOO)  (1) 

where: 

HTG2  =  height  increment  prediction  for  trees  with  d.b.h.  less  than  3  inches. 
HAB  =  constant  term  (intercept)  that  is  dependent  on  habitat  type. 
LOC    =  constant  term  (intercept)  that  is  dependent  on  location. 
ASP    =  stand  aspect  (degrees). 
SL      =  stand  slope  ratio  (%/100). 
HT      =  tree  height  (ft). 
CCF    =  crown  competition  factor. 
BAL    =  based  area  in  larger  trees  (ft^/acre). 
bpbg.bg  =  regression  coefficients  that  are  dependent  on  species  (see  tables  1 
and  2). 


Table  1.— Coefficients  for  the  small  tree  height  increment  model  (see  equation  1) 


 Species^  

Variable  WP  L  DF  GF  WH  C  LP  S  AF  PP  MH 


In(HT) 

(bi) 

0.4214 

0.2716 

0.3907 

0.3487 

0.3417 

0.2354 

0.5843 

0.2827 

0.3740 

0.4485 

0.2354 

CCF 

(bj) 

-  .0059 

-  .0065 

-  .0059 

-  .0039 

-  .0039 

-  .0039 

-  .0065 

-  .0039 

-  .0039 

- .0065 

-  .0039 

BAL/ 100 

(ba) 

-.3720 

-.4153 

-  .4004 

-  .2536 

-  .3469 

-.1201 

-.2417 

-.2530 

-  .2296 

-  .4730 

-  .2535 

Habitat 

1 

1.2554 

1.4058 

1 .2786 

.7835 

.8056 

.6807 

1.0190 

.8818 

.8521 

1.5165 

.6807 

class 

2 

1.3759 

1.5263 

1.3991 

.9040 

.9261 

.8012 

1.1392 

1.0023 

.9726 

1.6370 

.8012 

constants 

3 

1.1559 

1.2908 

.9531 

.7205 

1.0202 

.8953 

.9852 

.7533 

.5751 

1.2966 

.5215 

(HAB)2 

4 
5 
6 

1.4700 

1.6204 

1.0984 
1 .4932 

.9981 

.7202 
.8841 
1.2336 

.7085 
1.0667 

1.7311 

.8953 

Location 

1 

0.0 

class 

2 

-.0480 

[          same  for  all  species 

constants 

3 

-.2785  , 

(L0C)3 


'Species  codes  are  defined  in  appendix  A,  table  14. 
^Habitat  classes  are  defined  in  table  2. 
■'Location  classes: 

1  Clearwater  and  Nezperce  National  Forests 

2  St.  Joe  and  Coeur  d'Alene  National  Forests 

3  All  othier  Forests. 
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Table  2.  — index  Dy  species  to  hat3itat  class  constants  for  the  small  tree  height 
increment  model  (see  equation  ii 


Habitat   Spe^cies^' 


code^ 

WP 

L 

OF 

GF 

WH 

C 

LP 

S 

AF 

pp 

MH 

130 

3 

3 

i 

3 

5 

4 

1 

3 

170 

3 

3 

i 

3 

5 

4 

1 

3 

250 

3 

3 

3 

5 

4 

1 

3 

260 

3 

3 

A 

3 

5 

4 

3 

3 

280 

3 

3 

4 

3 

1 

4 

3 

3 

290 

3 

3 

i 

3 

1 

4 

1 

3 

310 

3 

3 

i 

3 

1 

4 

1 

3 

320 

3 

3 

1 

3 

5 

4 

3 

3 

330 

3 

3 

4 

3 

5 

4 

3 

3 

420 

3 

1 

4 

3 

5 

4 

3 

3 

470 

3 

1 

4 

3 

5 

4 

3 

3 

510 

3 

1 

1 

3 

1 

4 

1 

3 

520 

1 

1 

1 

1 

1 

1 

1 

1 

530 

2 

2 

2 

2 

2 

2 

2 

2 

2 

2 

2 

550 

4 

5 

4 

3 

3 

6 

4 

5 

4 

4 

570 

4 

4 

5 

4 

3 

3 

6 

4 

5 

4 

4 

610 

4 

4 

5 

4 

3 

3 

6 

4 

5 

4 

4 

620 

1 

1 

1 

1 

6 

'' 

1 

1 

1 

640 

3 

3 

4 

3 

4 

'' 

4 

3 

3 

660 

3 

3 

3 

3 

4 

4 

3 

3 

670 

3 

3 

4 

1 

3 

4 

3 

2 

680 

3 

3 

4 

1 

4 

4 

3 

3 

690 

3 

3 

4 

3 

5 

4 

3 

3 

710 

3 

3 

4 

3 

5 

4 

3 

3 

720 

3 

3 

4 

3 

4 

1 

3 

3 

730 

3 

3 

4 

3 

4 

1 

3 

3 

830 

3 

3 

3 

3 

4 

3 

3 

3 

850 

3 

3 

4 

3 

5 

3 

3 

3 

999 

3 

3 

4 

3 

5 

4 

3 

3 

Habitat  and  species  codes  are  defined  in  appendix  A.  t-ables  13  and  14 


In  addition,  the  strategy-  for  predicting  diameter  increments  for  small  trees 
was  changed.  In  version  4.0.  the  large-tree  d.b.h.  increment  model  is  applied  to 
small  trees  with  the  frequent  result  that  trees  with  unusual  proportions  are 
generated  (usually,  large  heights  with  small  d.b.h.  s).  The  new  model  is  applied 
only  to  trees  with  d.b.h.  less  than  3  inches,  and  it  predicts  d.b.h.  directly  from 
height  with  adjustments  for  stand  density  and  relative  size  (eq.  2).  The  adjust- 
ments were  developed  by  analyzing  prediction  errors  from  a  white  pine  spacing 
study  in  Deception  Creek  Experimental  Forest  near  Coeur  d  Alene.  ID. 

DBH  =  bo(HT-4.5)bi  -t- 

FAXl?  (0.01232-CCF  -  1.75)-RELh1-  (RELH  -  2.0)+  0.65  (2) 
L    36  J 

where: 

DBH  =  diameter  at  breast  height. 

AVH  =  average  height  of  the  40  trees  per  acre  with  the  largest  d.b.h.  s  (top 
height). 

RELH  =  relative  height  =  (HT-4.5)/(AVH-4.5)  (O  <  RELH  <  1.0). 
b^.b,  =  regression  coefficients  that  are  dependent  on  species  (see  table  3). 

Finally,  the  procedure  for  assigning  crown  ratio  (CR)  to  small  trees  was 
replaced.  The  new  model  is  based  on  the  logistic  equation  and  predictions  are 
dependent  on  tree  height,  tree  d.b.h.,  and  stand  basal  area  (eq.  3).  Parameters 
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Table  3.  — Coefficients  for  the  model  used  to  predict  d.b.h. 
(inches)  for  small  trees  (see  equation  2) 


Species 

bo 

bi 

Western  white  pine 

0.0781 

1 .1645 

Western  larch 

,0751 

1.1 176 

Douglas-fir 

.0828 

1.1713 

Grand  fir 

.1 155 

1.0688 

Western  hemlock 

.0729 

1.1988 

Western  redcedar 

.0730 

1.2343 

Lodgepole  pine 

.0988 

1.0807 

Engelmann  spruce 

.0658 

1.3817 

Subalpine  fir 

.0658 

1.3817 

Ponderosa  pine 

.2160 

1.0049 

Mountain  hemlock 

.0729 

1,1988 

Table  4.  — Coefficients  for  the  model  used  to  assign  crown  ratio  to 
small  trees  (see  equation  3) 


Species 

bo 

b. 

b2 

b3 

Western  white  pine 

-0.4432 

-0.4845 

0.0582 

0.00513 

Western  larch 

-  .8396 

-.1611 

.0416 

.00602 

Douglas-fir 

-.8912 

-.1808 

,0519 

.00454 

Grand  fir 

-  .6265 

-.0614 

,0236 

,00505 

Western  hemlock 

-.4955 

.0 

,0036 

,00456 

Western  redcedar 

.1185 

-.3931 

.0278 

,00626 

Lodgepole  pine 

-  .3247 

-.2011 

.0422 

,00436 

Engelmann  spruce 

-  .9201 

-.2245 

.0325 

,00620 

Subalpine  fir 

-  .8901 

-.1803 

.0223 

,00614 

Ponderosa  pine 

-.1756 

-.3385 

.0570 

.00692 

Mountain  hemlock 

-  .4955 

.0 

.0036 

.00456 

were  estimated  from  the  data  used  to  develop  the  small-tree  height  increment 
model.  For  most  species.  CR  increases  with  increasing  d.b.h.  and  decreases  with 
increasing  HT  or  BA.  Western  and  mountain  hemlock,  however,  are  insensitive 
to  changes  in  d.b.h. 


-     I  +  exp  (bo  +  bj-DBH  +  b,-HT  +  bg-BA) 
where: 

BA  =  stand  basal  area  (ft^/acre). 
bQ,bj.b2,b3  =  regression  coefficients  that  are  dependent  on  species  (see  table  4). 

If  the  inventory  does  not  include  some  or  all  crown  ratio  measurements,  this 
model  is  used  to  assign  CR  for  input  tree  records  that  are  less  than  3  inches 
d.b.h.  The  model  is  also  used  to  assign  a  new  CR  during  the  cycle  in  which  the 
tree  attains  a  3-inch  d.b.h.  Thereafter,  CR  change  is  predicted  as  in  version  4.0. 

Crown  ratio  does  not  influence  small  tree  growth  predictions  and,  therefore, 
no  periodic  change  is  simulated.  Missing  input  values  are  dubbed  to  provide 
consistency  in  program  output  and  to  facilitate  operation  of  Prognosis  Model 
extensions  such  as  COVER  and  the  Douglas-fir  tussock  moth  model. 

As  a  result  of  model  changes,  predictions  of  height  and  d.b.h.  increment  are 
sensitive  to  both  relative  tree  size  and  overall  stand  density  (fig.  1).  This  sys- 
tem of  models  has  demonstrated  reasonable  behavior  in  tests  to  date  although 
comparisons  with  data  from  larch  spacing  studies  (Milner  1985)  show  clear 
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VERSION  5.0 


INITIAL  HEIGHT,  FT 

Figure  1— Height  and  diameter  increment  predictions  for 
small  trees  are  sensitive  to  social  position  (as  indicated  by 
percentile  in  the  basal  area  distribution,  PCT)  and  to  overall 
stand  density  (crown  competition  factor,  CCF),  Curves 
shown  are  for  Douglas-fir  growing  on  a  western  redcedar/ 
CImtonia  habitat  type  at  3,700  ft  elevation  in  the  Clearwater 
National  Forest. 


DBH,  IN 

Figure  2.— The  most  notable  differences  in  large-tree  d.b.h. 
increment  predictions  are  for  subordinate  trees  growing  in 
dense  stands.  As  stand  density  increases,  and  as  percentile 
in  the  basal  area  distribution  decreases,  the  maximum  in 
the  increment  curve  shifts  downward  to  correspond  with 
larger  d.b.h.  Conditions  displayed  are  identical  to  those  dis- 
played in  Wykoff  and  others  (1982;  fig.  29,  p.  63). 


trends  in  prediction  errors  for  both  height  and  d.b.h.  increment  relative  to  ini- 
tial spacing.  Many  of  the  relationships  are  based  on  crude  analyses  or  limited 
data  or  both,  and  the  system  is  at  best  tenuous.  Work  is  currently  underway  to 
further  refine  the  prediction  of  small  tree  growth. 

Two  changes  were  made  in  the  version  5,0  large-tree  d.b.h.  increment  model 
(d.b.h.  >  3  inches).  The  ratio,  BAL/ln(DBH  -I-  1)  allows  the  magnitude  of  the 
relative  size  effect  to  vary  with  d.b.h.  In  addition,  the  CCF  coefficients  vary  by 
habitat  class.  As  a  result  of  these  changes,  the  model  is  less  sensitive  to  chang- 
ing d.b.h.  and  more  sensitive  to  stand  density  and  relative  size  (fig.  2). 


D.b.h.  Increment  for 
Large  Trees 
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In(dds)  =  HAB  +  LOG 

+  bi-SL-cos(ASP)  +  b2-SL-sin(ASP)  +  bg-SL  +  b4-SL2 
+  bg-EL  +  b6-EL2  +  b/lnlDBH) 

+  bg-CR  +  b9-CR2  +  bio-(BAL/100)  +  bu-BAL/lnfDBH  +  1) 

+  bi2-DBH2  +  bi3-(CCF/100)  (4) 

where: 

dds  =  10-year  change  in  squared  d.b.h. 
EL  =  stand  elevation  (in  hundreds  of  feet). 
CR  =  ratio  of  crown  length  to  total  tree  height, 
bj  through  b^3  =  regression  coefficients  that  are  dependent  on  species,  b^g  is 

dependent  on  location,  and  bj3  is  dependent  on  habitat  type 

(see  tables  5,  6,  7,  8,  and  9). 

Table  5.— Coefficients  for  the  large-tree  d.b.h.  increment  nnodel  (see  equation  4) 


Variables   Species 


WP 

W  r 

1 

U 1 

va  r 

vv  n 

1  P 

c 
o 

AP 

Mr 

pp 

IV1  n 

Habitat 

-| 

1.1558 

0.3834 

0.4778 

0.6676 

0.4526 

1 .6145 

0.7740 

-0.5884 

-  0  9639 

1.1623 

-1.6803 

class 

2 

1.0564 

.5129 

1523 

.6045 

1.3177 

.6783 

-  .2124 

-  .7242 

.7341 

- 1.521 1 

constants 

3 

.4538 

.2976 

.6445 

-  .7163 

-  .5731 

.5142 

(HAB)2 

4 
5 
6 

.7132 
.2684 

.3794 
.5434 

-  .5395 

-  .8222 
-1.2409 
-1.1075 

Location 

1 

.1692 

.2000 

.5036 

.4344 

.1067 

.5007 

.4374 

.2626 

.4206 

.2459 

.1252 

class 

2 

.0 

.0766 

.3492 

.2834 

.4436 

•1  ICC 

.1  ^bo 

.21 1 1 

^  C  0"7 

-  .100/ 

acne 

A  orto 

.4oUo 

constants 

3 

.0819 

.2196 

-.1483 

.0 

.3174 

.1481 

.0 

-.1300 

.4279 

.0 

{L0C)3 

4 
5 
6 

.3038 
.0 

.6181 
.0 

.2020 
.5776 
.0 

.0 

•0. 

.0 

.0 

SL»Cos(ASP) 

(bi) 

.0982 

-.2134 

-  .0456 

-.0122 

.0828 

-  .0662 

.0032 

-.1309 

-.1247 

-  .0998 

.1794 

SL«Sln(ASP) 

(bs) 

.0388 

.0343 

.0629 

-  .0460 

.1099 

.0553 

.1299 

-  .0604 

-  .0686 

.0119 

.1336 

SL 

(ba) 

-.1789 

.3352 

.7818 

1.1702 

.0497 

.1193 

.4655 

.6562 

.3007 

-  .0664 

.0763 

SL2 

(b4) 

.0 

-.7022 

-1.1238 

-1.5201 

.0 

.0 

-.5801 

-.9014 

-.6222 

-.4372 

.0 

EL 

(bs) 

.0352 

.0373 

.0259 

.0092 

.0286 

-.0018 

-.0048 

.0626 

.0631 

.0323 

.0852 

EL2 

(bg) 

-  .00047 

-  .00043 

-  .00038 

-.00012 

-  .00042 

- .00007 

-  .00006 

-.00071 

-  .00068 

- .00042 

-  .00094 

In(DBH) 

(b7) 

.5644 

.5414 

.5689 

.6881 

.6871 

.5870 

.8950 

.7304 

.8624 

.6610 

.8978 

CR 

(bg) 

1.0834 

1.0348 

2.0685 

1.9397 

1.6413 

1.2936 

1.8556 

1.5464 

.5204 

1.3162 

1.2840 

CR2 

(bg) 

.0 

.0751 

-.6236 

-.7826 

-.2724 

.0 

-  .3639 

-.2664 

.8624 

.0 

.0 

BAL/100 

(bio) 

.4211 

.4364 

.5020 

.4514 

.0 

.7460 

-  .0366 

.2564 

.0 

.0 

.0 

BAL/ln(DBH  +  1) 

(bi,) 

-2.0827 

-2.0326 

-2.1159 

-1.7681 

-  .8092 

-2.2838 

-.4333 

-1.1822 

-.5127 

-1.2588 

-.6611 

DBH2 

1 

-  .00044 

-  .00031 

-  .00025 

-  .00027 

-  .00022 

.0 

-.00126 

-.00013 

-  .00028 

-.00041 

-  .00048 

classes 

2 

.0 

-  .00057 

-  .00037 

-  .00009 

-  .00043 

-.00217 

-  .00029 

-  .00078 

-  .00044 

-.00031 

3 

-  .00050 

-  .00064 

-.00189 

-.00043 

-.00014 

4 

-  .00057 

-.00087 

CCF 

1 

-.0243 

-.1014 

-.0905 

-.0962 

.0 

-  .0505 

-.0558 

-.0155 

-.0160 

-.1042 

-.1074 

classes 

2 

-.2489 

-.1479 

-.1188 

-.1954 

-.1536 

-.1492 

-  .3839 

-.0448 

-  .8881 

3 

-.0108 

-  .0544 

-  .0553 

-.0512 

-  .0940 

-  .4064 

-.0537 

-.0739 

-  .2594 

4 

-.0218 

-.1140 

-.1516 

-.1473 

''Species  codes  are  defined  in  appendix  A,  table  14. 
^Habitat  classes  are  defined  in  table  6. 
-'Location  classes  are  defined  in  table  7. 
■•DBH-squared  classes  are  defined  in  table  8. 
^CCF  classes  are  defined  in  table  9. 
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Table  6.  — Index  by  species  for  the  habitat  constants  in  the  large-tree  d.b.h. 
increment  model  (see  equation  4) 


Habitat   Species^ 


codo^ 

WP 

L 

DF 

GF 

WH 

c 

LP 

s 

AF 

PP 

MH 

130 

2 

5 

3 

2 

2 

5 

4 

6 

1 

2 

170 

2 

5 

3 

2 

2 

5 

4 

6 

1 

2 

250 

2 

5 

3 

2 

2 

5 

4 

5 

2 

2 

260 

2 

5 

3 

2 

2 

5 

4 

5 

3 

2 

280 

2 

5 

3 

2 

2 

1 

4 

6 

3 

2 

290 

2 

5 

3 

2 

2 

2 

4 

5 

2 

2 

310 

2 

5 

3 

2 

2 

1 

4 

5 

2 

2 

320 

2 

5 

1 

2 

2 

5 

4 

6 

3 

2 

330 

2 

5 

3 

2 

2 

5 

4 

5 

3 

2 

420 

2 

1 

3 

2 

2 

5 

4 

6 

3 

2 

470 

2 

1 

3 

2 

2 

5 

4 

6 

3 

2 

510 

2 

2 

2 

2 

2 

1 

5 

2 

2 

520 

1 

1 

1 

2 

2 

1 

1 

2 

2 

530 

1 

2 

2 

2 

3 

4 

2 

2 

2 

550 

1 

3 

2 

1 

3 

2 

3 

1 

2 

570 

1 

3 

2 

2 

3 

4 

4 

1 

2 

610 

1 

3 

2 

2 

3 

2 

3 

1 

2 

620 

1 

2 

2 

2 

3 

1 

1 

3 

2 

640 

2 

5 

3 

2 

2 

4 

4 

6 

3 

2 

660 

2 

2 

2 

2 

2 

4 

4 

6 

3 

2 

670 

1 

1 

3 

I 

1 

3 

4 

5 

3 

1 

680 

1 

1 

3 

2 

2 

4 

4 

6 

3 

2 

690 

2 

1 

3 

2 

2 

5 

4 

6 

3 

2 

710 

2 

5 

3 

1 

2 

5 

4 

6 

3 

2 

720 

2 

4 

3 

2 

2 

4 

4 

1 

3 

2 

730 

2 

4 

3 

2 

2 

4 

4 

1 

3 

2 

830 

2 

5 

2 

2 

2 

4 

3 

5 

3 

2 

850 

2 

5 

3 

2 

2 

5 

4 

5 

3 

2 

9992 

2 

5 

3 

2 

2 

5 

4 

6 

3 

2 

^Habitat  and  species  codes  are  defined  in  appendix  A,  tables  13  and  14. 
^ypes  grouped  with  999  were  Included  in  the  overall  mean  for  the  species. 


Table  7.— Index  by  species  for  the  location  constants  in  the  large-tree  d.b.h, 
increment  model  (see  equation  4) 


National   Species^ 


Forest 

WP 

L 

DF 

GF 

WH 

C 

LP 

S 

AF 

PP 

MH 

Bitterroot 

2 

1 

5 

6 

3 

4 

4 

3 

4 

1 

3 

Clearwater 

2 

1 

1 

1 

3 

1 

1 

1 

1 

2 

1 

Coeur  d'Alene 

2 

2 

2 

2 

1 

1 

1 

1 

2 

2 

1 

Colville 

2 

3 

3 

2 

3 

2 

2 

3 

2 

1 

3 

Flathead 

2 

3 

3 

3 

3 

2 

4 

2 

3 

4 

3 

Kaniksu 

2 

2 

2 

2 

3 

3 

3 

3 

3 

3 

3 

Kootenai 

2 

5 

3 

4 

3 

4 

3 

3 

4 

1 

3 

Lolo 

2 

5 

5 

6 

3 

2 

4 

3 

4 

4 

1 

Nezperce 

2 

4 

1 

2 

3 

1 

2 

1 

2 

3 

3 

St.  Joe 

1 

1 

4 

5 

2 

1 

2 

1 

1 

2 

2 

^Species  codes  are  defined  in  appendix  A.  table  14. 
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Table  8.— Index  by  species  for  the  location-dependent  DBH-squared  coefficients  in 
the  large-tree  d.b.h.  increment  model  (see  equation  4) 


National  Species^ 


Forest 

WP 

L 

DF 

GF 

WH 

C 

LP 

S 

AF 

PP 

MH 

Bitterroot 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

Clearwater 

2 

2 

1 

2 

2 

2 

2 

2 

Coeur  d'Alene 

2 

1 

2 

1 

1 

2 

1 

1 

2 

1 

Colville 

2 

2 

1 

1 

1 

1 

2 

2 

Flathead 

1 

3 

2 

1 

1 

1 

1 

3 

Kaniksu 

2 

1 

2 

1 

2 

3 

1 

3 

Kootenai 

1 

4 

3 

2 

3 

2 

2 

2 

Lolo 

1 

1 

1 

1 

1 

1 

1 

1 

2 

Nezperce 

1 

1 

2 

2 

4 

1 

1 

1 

1 

St.  Joe 

2 

2 

4 

1 

2 

1 

1 

1 

2 

2 

2 

^Species  codes  are  defined  in  appendix  A,  table  14. 


Table  9.— Index  by  species  for  the  habitat-dependent  CCF  coefficients  in  the  large- 
tree  d.b.h.  increment  model  (see  equation  4) 


Habitat   Species^ 


code^ 

WP 

L 

DF 

GF 

WH 

C 

LP 

S 

AF 

PP 

MH 

130 

3 

3 

4 

3 

3 

4 

4 

3 

2 

O 

o 

A 

O 

O 

O 

A 

^ 

A 
H 

O 

o 

c. 

250 

3 

3 

4 

3 

3 

1 

4 

3 

3 

260 

3 

3 

4 

3 

3 

4 

4 

3 

1 

280 

3 

3 

4 

3 

3 

3 

4 

3 

4 

290 

3 

3 

4 

3 

3 

4 

4 

3 

3 

310 

3 

3 

1 

3 

3 

2 

4 

3 

3 

320 

3 

3 

2 

3 

3 

1 

4 

3 

1 

330 

3 

3 

4 

3 

3 

4 

4 

3 

2 

420 

3 

3 

4 

3 

3 

1 

4 

3 

4 

470 

3 

3 

4 

3 

3 

4 

4 

3 

4 

510 

3 

3 

2 

3 

3 

4 

3 

3 

1 

520 

1 

1 

1 

1 

3 

4 

3 

1 

4 

530 

3 

3 

4 

3 

1 

2 

1 

1 

4 

550 

3 

3 

4 

3 

2 

2 

1 

1 

4 

570 

1 

3 

3 

3 

3 

2 

3 

1 

1 

610 

3 

3 

4 

3 

3 

2 

3 

1 

4 

620 

3 

2 

4 

3 

2 

2 

1 

1 

4 

640 

3 

3 

4 

3 

3 

4 

4 

3 

4 

660 

3 

1 

1 

3 

3 

1 

2 

1 

4 

670 

2 

3 

3 

3 

1 

4 

4 

1 

4 

680 

2 

3 

2 

3 

3 

1 

1 

2 

4 

690 

3 

1 

4 

2 

3 

4 

4 

2 

4 

710 

3 

1 

4 

2 

3 

4 

4 

1 

4 

720 

3 

3 

4 

3 

3 

1 

1 

1 

4 

730 

3 

3 

4 

3 

3 

1 

1 

1 

4 

830 

3 

3 

2 

3 

3 

4 

4 

1 

4 

850 

3 

3 

4 

3 

3 

4 

4 

3 

4 

999 

3 

3 

4 

3 

3 

4 

4 

3 

4 

'Habitat  and  species  codes  are  defined  in  appendix  A,  tables  13  and  14. 
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Individual  Tree  The  mortality  model  now  incorporates  explicit  terms  that  reflect  stand  den- 

Mortality  Models  sity  and  tree  vigor.  In  version  4.0,  mortality  rate  was  predicted  from  a  regres- 

sion model,  with  DBH  and  DBH-  used  as  independent  variables.  Then,  various 
correction  factors  were  applied  to  reflect  the  effect  of  stand  density.  The  correc- 
tion factors  were  derived  from  yield  tables  (Haig  1932),  published  data  on  car- 
rying capacities  of  various  habitat  types  (Pfister  and  others  1977;  Daubenmire 
and  Daubenmire  1968),  analysis  of  Northern  Region  timber  management  plan- 
ning inventories,  and  "best  guesses."  The  stand  level  models  were  required 
because  the  data  used  to  develop  the  individual  tree  models  contained  limited 
information  on  stand  density  and  relative  size,  and  these  effects  were  therefore 
absent  from  the  fined  model. 

The  mortality  model  used  in  version  5.0  uses  rate  estimates  (R^  and  Rj,)  that 
are  predicted  from  two  independent  equations.  R^^  (eq.  5)  is  based  primarily  on 
the  anadysis  of  permanent  sample  plots  that  are  maintained  by  the  Intermoun- 
tain  Station  (Hamilton  in  preparation  a).  Predictions  are  dependent  on  habitat 
type,  species,  d.b.h.,  d.b.h.  increment,  estimated  potential  d.b.h.  increment, 
stand  basal  area,  and  relative  diameter  (d.b.h. /mean  stand  d.b.h.). 

R ,  =  11  +  exp[bo  +  b,-\  DBH  +  b^-v^X^  h^-g  +  b^-RDBH  + 

(b,  +  b,-g)-DBH~'  ]  I  -1  (5) 
where: 

Rg  =  estimated  annual  mortality  rate. 

g  =  periodic  annual  d.b.h.  increment  for  previous  growth  period  adjusted 
for  differences  in  potential  annual  d.b.h.  increment  indexed  by  habi- 
tat type  and  National  Forest. 
RDBH  =  the  ratio  of  tree  d.b.h.  to  the  arithrnetic  mean  stand  d.b.h. 

b„  =  species  dependent  constant  (see  table  10). 

b,  =  0.2223 

bg  =  -0.0460 

bg  =  10.0810 

b.,  =  0.2463 

b-  =  -0.5544 

bg  =  6.0713. 

Rg  is  multiplied  by  a  factor  based  on  Reineke's  (1933)  Stand  Density  Index  that 
accounts  for  expected  differences  in  mortality  rates  on  different  habitat  types 
and  National  Forests  (Hamilton  in  preparation  b). 

Table  10.  — Species-dependent 

constants  for  the  mortality 
model  (see  equation  9) 


Species  Constant  (bg) 


Western  white  pine 

0.0 

Western  larch 

-.1760 

Douglas-fir 

.3179 

Grand  fir 

.3179 

Western  hemlock 

.6077 

Western  redcedar 

1.5798 

Lodgepole  pine 

-.1206 

Engelmann  spruce 

.9402 

Subalpine  fir 

.2118 

Ponderosa  pine 

.2118 

f\/lountain  hemlock 

.0 
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DOUGLAS-FIR:   VERSION  l4 .  1 


SYMBOL   PRINTED  REPRESENTS  THE  NUMBER  OF 
DECADES   PROJECTED  PRIOR  TO  THE   PERIOD  FOR 
WHICH   THE  RATE   ESTIMATE  WAS  MADE. 

6  +PREDICTED  ANNUAL 

IMORTALITY  RATE  (PERCENT)  PROJECTED       STAND  BASAL 

DECADES         AREA  (SQFT) 

0  7 

1  39 
3  128 
5  214 
7  258 
9  295 


4 


3 


00 
0000 
000 
100 

1111 
1 1 1 

1 1 

1 

1 


3 

3333  5  7777777  7777  77     7  7     7777  5  577  7 

33  5  5     55  5    9      55    95        99  59999  9        9  9 
3  333 


0  +  +  ---+   +   *   +  

0  3  6  9  12  15  18 

DIAMETER  BREAST  HEIGHT  (IN. 


DOUGLAS-FIR:   VERSION  5.1 

SYMBOL  PRINTED  REPRESENTS  THE  NUMBER  OF 
DECADES  PROJECTED  PRIOR  TO  THE   PERIOD  FOR 
WHICH  THE  RATE  ESTIMATE  WAS  MADE. 

PREDICTED  ANNUAL 

MORTALITY  RATE   (PERCENT)  PROJECTED       STAND  BASAL 

DECADES         AREA  (SQFT) 


0  7 

1  65 
3  205 
5  277 
7  300 
9  285 


0  3 

3  77 

5  7 

9  5 

7  7 
7  7 
0  77  7 

0  5         5  55  99  9 

0  7  9  9 

0  9       5  7 

0  1  3  5  7  9 

0  3     5  7  7  9 

0  359  9  9     9-999  7 

0     1  5  57  5  5  9      77  9 

0  3  5     755       7       7  9 

0  1  3         3  7  9 

0  3  35  9     9  99 

0  3  3  3  555      97  9        5     57  7     5  7  9 

0  3  3  555955579  7799 

01  33         333       33         5777  55  555 

0     1  3  3  3       3  3  7 

0       1  3  33       333         333     33  77  7     97       7  9  99 

00       1  73  3         93  7     7     9  9  9  7  7 

0  111 

00  1    nil  1 

no       111111  1 

11111111  1 
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DIAMETER  BREAST  HEIGHT  (IN.) 

Figure  3.  — Predicted  annual  mortality  rate  versus  d.b.h.  for  Douglas-fir.  The  numbers 
printed  represent  decades  in  a  long-term  projection  (1  is  the  average  rate  for  years 
11-20  in  the  projection;  9  is  the  average  rate  for  years  90-100  in  the  projection).  Ini- 
tial stand  and  site  conditions  are  the  same  as  for  figure  1.  Figure  3a  displays  rates 
projected  by  version  4.1;  figure  3b  displays  rates  projected  by  version  5.1. 
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For  purposes  of  predicting  mortality  rate,  d.b.h.  and  d.b.h.  increment  are 
bounded.  If  d.b.h.  is  less  than  0.5  inch,  it  is  set  equal  to  0.5  inch:  if  d.b.h.  is 
less  than  1  inch  and  d.b.h.  increment  is  less  than  0.05  inch,  d.b.h.  increment  is 
set  to  0.05  inch. 

The  second  part  of  the  mortality  rate  estimate,  Rj,,  is  dependent  on  the  prox- 
imity of  stand  basal  area  to  the  assumed  maximum  for  a  site  (BAMAX)  and  on 
the  estimated  rate  of  basal  area  increment.  As  stand  basal  area  approaches 
BAMAX.  R(j  approaches  1.  The  calculation  of  Rj,  is  described  by  Wykoff  and 
others  (1982).  As  in  version  4.0,  the  value  of  BAMAX  is  indexed  to  habitat 
type,  but  the  default  value  may  be  overridden  by  the  user. 

The  mortality  rate  actually  applied  to  a  tree  is  a  weighted  average  of  R^  and 
R^,.  The  weights  applied  to  the  respective  estimates  are  also  dependent  on  the 
proximity  of  stand  basal  area  to  BAMAX  (eq.  6). 

R,  =  W-R^  +  (l-W)-R^  (6) 
where: 

—  annual  mortality  rate  applied  to  tree  t. 
W  =  BABAMAX. 

In  general,  the  new  mortality  model  is  more  sensitive  to  stand  density  and  to 
the  distribution  of  size  classes  (fig.  3)  and  better  reflects  the  relationships 
between  stand  dynamics  and  management. 

Random  Effects  In  application,  random  errors  are  drawn  and  attached  to  predictions  of  large- 

tree  d.b.h.  increment,  small-tree  height  increment,  and  small-tree  crown  ratio. 
For  all  of  these  variables,  the  error  is  assumed  to  be  Normally  distributed  when 
the  model  is  transformed  so  that  it  is  linear  in  the  parameters.  The  random 
error  is  drawn  from  a  Normal  distribution  on  the  transformed  scale.  Random 
effects  are  propagated  in  other  models  by  including  d.b.h.,  crown  ratio,  or 
height  increment  or  all  three  as  independent  variables. 

When  nonstockable  plots  are  encountered  in  a  stand  inventory,  their  ratio  to 
total  plots  is  taken  to  represent  the  proportion  of  total  stand  eirea  that  is  unin- 
habitable by  trees.  In  version  4.0,  only  the  value  of  crown  competition  factor  is 
adjusted  to  reflect  stockability,  resulting  in  an  inconsistent  representation  of 
density  effects  in  the  increment  models.  In  version  5.0,  the  stand  area  used  to 
compute  the  trees  per  acre  represented  by  each  sample  tree  is  based  only  on 
stockable  plots.  As  a  result,  all  stand  and  tree  attributes  carried  by  the  pro- 
gram represent  stockable  area.  As  described  in  the  section  titled  "Changes  in 
Program  Output.  "  the  program  output  now  reports  stand  attributes  on  the 
basis  of  both  total  and  stockable  area. 

In  version  5.0.  there  is  a  new  procedure  for  using  input  gn'owth  data  to  adjust 
the  large-tree  d.b.h.  increment  and  small-tree  height  increment  models.  The  new 
procedure  is  designed  to  eliminate  a  bias  that  occurred  in  version  4.0  when  all 
sample  trees  were  measured  for  increment.  Predictions  are  made  for  each 
growth-sample  tree  and  then  observed  increment  (g^)  is  regressed  against 
predicted  increment  (gp): 

go  =  a  +  b-gp 

A  prediction  is  made  for  every  tree  to  obtain  a  mean  predicted  value  for  the 
population  (Gp).^The  regression  model  is  used  to  estimate  mean  increment  for 
the  population  (G^)  from  the  difference  between  the  mean  predictions  for  the 
population  and  for  the  growth  sample  (^p): 


Interpretation  of 
Input  Stockability 
Data 


Calibration  of 
Growth  Models 
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Attenuation  of 
Model  Calibration 


Assigning  Initial 
Values  for  Missing 
Heights 


Computing  Scribner 
Board  Foot  Volume 
with  Variable 
Merchantability 
Standards 


G„  =       +  b.(Gp  -  gp) 

The  ratio,  GJG^  is  then  used  as  a  multiplicative  adjustment  to  the  growth 
model. 

As  in  version  4.0,  the  model  for  each  species  is  cahbrated  independently.  In 
version  5.0,  however,  there  must  be  five  (instead  of  two)  increment  observations 
before  the  model  for  a  species  will  be  adjusted.  The  trees  on  which  d.b.h.  incre- 
ment is  sampled  must  be  greater  than  3  inches  d.b.h.  at  the  start  of  the  growth 
period. 

The  model  adjustment  that  is  based  on  input  increment  data  is  attenuated 
over  time.  For  both  the  large-tree  d.b.h.  increment  model  and  the  small-tree 
height  increment  model,  the  asymptote  for  attenuation  is  defined  as  the  aver- 
age of  1.0  and  the  original  adjustment  to  the  d.b.h.  increment  model.  An 
exponential  decay  function  is  then  used  to  approach  the  asymptote  over  time;  a 
value  midway  between  the  original  adjustment  and  the  asymptote  will  be 
attained  25  years  after  the  start  of  the  projection. 

When  the  height  for  a  tree  is  missing  (as  indicated  by  a  blank  or  zero  input 
value),  a  height  is  assigned  from  a  height-diameter  function: 

HT  =  exp[Co  +  Ci-1/(DBH  +  1)]  +  4.5 

The  model  parameters  are  given  in  the  User's  Guide  to  the  Stand  Prognosis 
Model  (Wykoff  and  others  1982,  p.  52). 

Version  4.0  used  a  regression  procedure  to  reestimate  Cq  and  Cj  when  there 
were  four  or  more  height  observations  for  a  species.  In  version  5.0,  however,  Cj 
is  held  constant  and  the  height  sample  is  used  to  adjust  Cq.  Thus  the  level  of 
the  curve  shifts  but  shape  is  retained.  By  limiting  the  calibration  to  the  inter- 
cept term,  the  procedure  is  less  vulnerable  to  a  poor  choice  of  height  sample 
trees.  If  for  instance,  the  height  sample  contained  only  large  trees  with  uniform 
heights,  the  estimate  for  Cj  could  be  positive,  resulting  in  unrealistically  tall 
height  estimates  for  trees  with  small  d.b.h. 's. 

Volume  calculations  have  been  modified  in  version  5.0  to  permit  specification 
of  variable  top  diameter  for  computing  board  foot  volume.  The  modification 
was  developed  by  J.  E.  Brickell,  mensurationist  with  the  USD  A  Forest  Service, 
Northern  Region  Timber  Management  staff,  Missoula,  MT. 

Version  4.0  used  equations  developed  by  P.  D.  Kemp  (see  Wykoff  and  others 
1982,  p.  82),  with  fixed  merchantability  standards  (9-inch  minimum  d.b.h.,  1-foot 
stump  height,  and  8-inch  minimum  top  diameter).  An  additional  set  of  volume 
equations  (Allen  and  others  1974)  predicts  board  foot  volume  assuming  a  2-inch 
minimum  d.b.h.  and  a  2-inch  minimum  top  diameter.  These  two  sets  of  equa- 
tions can  predict  volume  to  a  variable  top  diameter  (D^)  if  the  tree  is  considered 
as  four  segments  (fig.  4).  Each  segment  is  assumed  to  be  a  tree;  d.b.h.  of  each 
assumed  tree  is  estimated  from  the  known  basal  diameter  of  the  segment  and 
height  is  estimated  with  a  height-diameter  function. 

The  first  tree  segment  (Sj)  extends  from  a  1-foot  stump  height  to  a  height 
corresponding  to  a  that  is  greater  than  8  inches.  The  second  segment  (Sg) 
extends  from  the  top  of  S^  to  the  height  at  which  is  equal  to  8  inches.  The 
third  segment  (S3)  extends  from  the  top  of  Sg  to  a  height  at  which  is  less 
than  8  inches  but  greater  than  2  inches.  The  final  segment  (S4)  extends  from 
the  top  of  S3  to  the  height  at  which      is  equal  to  2  inches. 
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Table  11.  — Map  of  species  onto  volume  equations  for  computing 
Scribner  board  foot  volume  from  squared  d.b.h.  and 
IN  height  ^-  (D^H).  Allen's  equations  are  used  when  D^H  is 

less  than  the  transition  value.  Kemp  s  equations  are  used 
when  D^H  is  greater  than  the  transition  value 


Equation  used  Transitional  value 

Species  name  (Allen  and  others  1974)       of  D^H  (in^  x  ft) 


8  IN 


Western  white  pine 

Lodgepole  pine 

16.917 

Western  larch 

Western  larch 

19.232 

Douglas-fir 

Douglas-fir 

15.764 

Grand  fir 

Grand  fir 

10.017 

Western  hemlock 

Grand  fir 

14,884 

Western  redcedar 

Western  larch 

29.922 

Lodgepole  pme 

Lodgepole  pine 

4.568 

Engelmann  spruce 

Grand  fir 

6,025 

Subalpine  fir 

Grand  fir 

19.426 

Ponderosa  pme 

Douglas-fir 

14.034 

Mountain  hemlock 

Grand  fir 

14,884 

8  IN 


Figure  4.  — Diagram  of 
a  tree  showing  the 
four  tree  segments 
(S,.  Sj.  S3,  and  S^) 
that  are  used  to  com- 
pute Scribner  board 
foot  volume  to  varia- 
ble merchantability 
standards  using  two 
independent 
equations. 


Allen's  and  Kemp's  equations  are  both  dependent  on  the  product  of  squared 
diameter  and  height  (D^H).  For  each  species,  the  equations  intersect:  the  point 
of  intersection  (T)  corresponds  to  a  positive  value  of  D^H.  For  D^H  less  than  T, 
Alien  s  equations  are  used  to  predict  volume  to  a  2-inch  top  (S,  +  S2  +  S.3  + 
S,).  If  a      greater  than  2  inches  is  specified,  S,,  is  computed  with  Allen's  equa- 
tion and  deducted  from  the  total. 

For  D^H  greater  than  or  equal  to  T,  Kemp  s  equation  is  used  to  predict  vol- 
ume to  an  8-inch  top  (S,  +  S.^).  If  a      greater  than  8  inches  is  specified,  Kemp's 
equation  is  used  again  to  predict  the  volume  of  S^,  which  is  subtracted  from 
the  total.  If  a      of  less  than  8  inches  is  specified.  (S3  +  S^)  and  S4  are  esti- 
mated with  Allen's  equation  in  order  to  find  top  volume  (S3),  which  is  added  to 
the  total. 

Allen's  equations  were  developed  only  for  lodgepole  pine,  Douglas-fir,  western 
larch,  and  grand  fir.  Other  species  were  mapped  onto  one  of  these  four  equa- 
tions for  use  in  the  Prognosis  Model.  The  mapping  and  the  vEilue  of  D^H  at 
which  Allen's  and  Kemp's  equations  give  identical  results  are  given  in  table  11. 


13 


Tree  Record 
Compression 


Seeding  the  Random 
Number  Generator 


Generating  a  List  of 
Projected  Tree 
Records 


Version  5.0  contains  a  procedure  for  classifying  records  into  groups  that  are 
similar  with  regard  to  variables  that  strongly  influence  growth  predictions.  The 
purpose  of  this  compression  algorithm  is  threefold:  (1)  to  provide  space  for  new 
records  when  the  regeneration  establishment  model  is  called;  (2)  to  eUminate 
records  which,  due  to  thinning  or  mortality,  no  longer  represent  a  significant 
contribution  to  the  stand  as  a  whole;  and  (3)  to  reduce  the  number  of  tree 
records  so  that  yield  projections  are  less  costly  from  a  computational  stand- 
point, but  retain  much  of  the  original  variation  in  tree  characteristics. 

A  two-step  algorithm  is  used  to  compress  the  tree  records  (Stage  and  others 
in  preparation).  First,  a  classification  function,  which  is  a  weighted  linear  com- 
bination of  selected  tree  attributes  (including  species,  d.b.h.,  predicted  d.b.h. 
increment,  height,  crown  ratio,  and  basal  area  percentile)  is  used  to  find  the 
largest  gaps  in  a  one-dimensional  space.  These  gaps  are  used  to  identify 
clusters  of  similar  records.  In  the  second  step,  the  larger  of  the  identified 
clusters  (based  on  the  maximum  difference  in  the  value  of  the  classification 
function)  are  spht.  Following  the  second  step,  the  attributes  of  all  of  the 
records  in  a  cluster  are  averaged  to  make  a  single  record.  Unless  the  tEirget  is 
changed  by  the  user  (see  discussion  of  the  COMPRESS  keyword)  the  tree  list 
will  be  compressed  to  150  records. 

CHANGES  IN  FEATURES  THAT  WERE  PART  OF  VERSION  4.0 

In  version  4.0,  the  random  number  generator  was  initialized  only  at  the 
beginning  of  the  first  of  several  projections  in  a  runstream.  As  a  result,  chang- 
ing the  order  of  projections  within  a  runstream  changed  the  results  of  the 
projections.  In  version  5.0,  the  random  number  generator  is  initialized  before 
each  projection  so  that  repeatability  of  projections  is  not  confounded  with  their 
order  in  the  runstream.  The  sequence  of  random  numbers  can  be  changed  by 
using  the  RANNSEED  record  to  input  a  new  seed.  In  version  4.0  there  are 
three  seeds;  the  version  5.0  generator  requires  only  one. 

RANNSEED  field  1:    New  seed  for  random  number  generator.  The  seed 

should  be  a  positive,  odd  integer  value.  If  otherwise, 
fractions  will  be  truncated  and/or  1  will  be  added; 
default  =  55329. 

Once  the  RANNSEED  record  has  been  used  to  reseed  the  random  number 
generator,  the  input  seed  will  be  used  to  initialize  the  generator  in  subsequent 
projections  in  the  same  runstream. 

As  in  version  4.0,  a  complete  list  of  tree  records  can  be  generated  for  any  or 
all  cycles.  In  addition,  the  user  can  now  control  the  dataset  reference  number 
for  the  output  file  and/or  suppress  headings.  When  the  no-heading  option  is 
selected,  a  single  identification  record  is  printed  at  the  beginning  of  the  Ust  for 
each  cycle.  The  record  contains  a  flag  (the  number  —999  in  columns  1-4),  the 
number  of  records  printed  (columns  5-8),  and  the  cycle  number  (columns  12-13). 

TREELIST     field  1:    The  cycle  in  which  a  complete  Ust  of  trees  is  to  be 

printed.  The  hst  is  printed  at  the  end  of  the  cycle  and 
the  records  are  updated  to  include  growth  for  the 
period.  If  blank,  a  list  will  be  generated  at  the  begin- 
ning of  the  projection  and  at  the  end  of  each  cycle, 
field  2:    Dataset  reference  number  for  the  output  tree  record  file; 
default  =  3. 
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field  3:    Heading  suppression  indicator.  Any  numeric  entry  in 
this  field  will  prevent  headings  from  being  printed; 
default  =  print  headings. 

It  is  possible  to  produce  more  than  one  tree  list  in  any  cycle.  For  example,  if 
both  a  printed  list  and  a  machine-readable  list  are  desired,  a  second  TREELIST 
record  can  be  inserted  in  the  runstream  to  provide  a  dataset  reference  number 
for  the  machine-readable  file: 


TREEL  i  ST 
TREEL I  ST 


2  .  0 
2  0 


V) 

99.0 


9  9  9.0 


In  this  case,  the  files  will  be  written  at  the  end  of  cycle  2  and  the  second  file 
will  be  directed  to  unit  99  with  headings  suppressed. 


Additional 
Flexibility  in 
Field  Sampling 
Procedures 


Tree  Records  Can 
Be  Completely 
Removed  in 
Thinning 


The  interpretation  of  the  DESIGN  record  has  been  changed  so  that  it  is  now 
possible  to  sample  with  two  sizes  of  fixed-area  plots  where  choice  of  plot  size  is 
dependent  on  tree  d.b.h.  For  example,  a  1/20-acre  plot  could  be  used  to  sample 
large  trees  with  a  1/300-acre  plot  used  to  sample  seedlings  and  saplings.  A 
negative  value  entered  in  field  1  of  the  DESIGN  record  is  interpreted  as  the 
inverse  of  the  plot  area  from  which  large  trees  were  sampled.  Thus,  the  above 
example  could  be  specified  by: 

DF.S  IGN  -20  .0  300  .0  5  .0  16  Ift  16 

where  5  inches  is  used  as  the  breakpoint  (BRK)  for  determination  of  plot  size. 

DESIGN        field  1:    Positive  value  is  interpreted  as  basal  area  factor  for 

horizontal  angle  gauge;  negative  value  is  interpreted  as 
the  inverse  of  large-tree  fixed  plot  area;  default  =  40 
(square  feet/acre/tree), 
field  2:    Inverse  of  small-tree  fixed  plot  area;  default  =  300 
(acre" '). 

field  3:    BRK;  default  =  5  (inches). 

field  4:    Number  of  plots  in  the  stand.  If  blank,  or  zero,  the 

number  of  plots  in  the  stand  is  determined  by  counting 
the  numbers  of  unique  plot  identification  codes  on  the 
tree  records. 

field  5:    Number  of  nonstockable  plots  in  the  stand.  These 
include  plots  falling  on  rock  outcroppings,  roads, 
streams,  etc.  If  blank,  count  nonstockable  plots  on  tree 
records  (IMC  =  8;  see  discussion  of  tree  records  in  ver- 
sion 4.0  User's  Guide  [Wykoff  and  others  1982]). 

field  6:    Sampling  weight  for  stand.  This  weight  does  not  affect 
the  projection  but  is  used  in  programs  that  aggregate 
many  projections  to  produce  a  composite  yield  table; 
default  =  number  of  plots. 

In  version  4.0,  no  more  than  99  percent  of  any  tree  record  could  be  removed 
in  any  one  thinning.  This  constraint  has  been  removed.  The  CUTEFF  record  in 
version  5.0  can  specify  a  cutting  efficiency  parameter  ranging  from  0.01  to  1.0. 
If  1.0  is  specified,  and  a  tree  is  removed,  the  tree  record  will  be  deleted  from 
the  list.  The  maximum  of  1.0  also  applies  to  the  cutting  efficiency  parameter 
entered  on  a  specific  thinning  request  (see  next  page). 
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D.b.h.  Limits  Added 
to  Thinning  Options 


Access  to  the 
COVER  Extension 


CUTEFF        field  1:    Proportion  of  the  trees  per  acre  represented  by  a  record 

that  is  eliminated  if  a  tree  is  designated  for  removal  in 
any  thinning.  The  value  of  this  parameter  must  fall 
between  0.01  and  1.0  or  the  keyword  will  be  ignored; 
default  =  0.98. 

The  thinning  keywords  that  are  used  to  specify  a  stand  density  target  can 
now  be  applied  within  user-selected  d.b.h.  limits.  Only  trees  with  d.b.h. 's 
between  the  selected  limits  will  be  removed  in  a  thinning.  The  minimum  and 
maximum  d.b.h.  values  are  entered  respectively  on  fields  4  and  5  of  the  key- 
word records. 

THINBTA  —    Thin  from  below  to  a  trees-per-acre  target. 
THINATA  —    Thin  from  above  to  a  trees-per-acre  target. 
THINBBA  —   Thin  from  below  to  a  basal  area-per-acre  target. 
THINABA  —    Thin  from  above  to  a  basal  area-per-acre  target. 

field  1:    Year  in  which  thinning  is  requested;  if  blank,  schedule 

at  start  of  projection, 
field  2:   The  desired  residual  stand  density  (trees  per  acre  or 

basal  Eirea  per  acre,  depending  on  the  keyword).  If  a 

residual  density  is  not  specified,  the  thinning  request 

will  be  ignored. 

field  3:    The  cutting  efficiency  parameter  to  be  used  only  with 
this  thinning  request.  If  blank,  or  outside  the  range 
0.01  -  1.0,  use  the  value  specified  on  the  CUTEFF 
record. 

field  4:    Smallest  d.b.h.  to  be  considered  for  removal;  default  = 
0.0. 

field  5:    Largest  d.b.h.  to  be  considered  for  removal;  default  = 
999.0. 

As  an  example,  a  stand  is  to  be  thinned  from  above  to  a  residual  density  of 
100  trees  per  acre.  In  addition,  no  trees  that  are  less  than  16  inches  d.b.h.  are 
to  be  removed.  This  prescription  can  be  simulated  by: 

TH  I  NATA  1984.0  100.0  0.5  16.0  \b 

The  d.b.h.  limits  (16.0  and  999.0  in  this  example)  and  the  cutting  efficiency 
parameter  (0.5  in  this  example)  may  constrain  the  actual  removal  to  less  than 
the  specified  target.  If,  in  this  example,  there  are  200  trees  per  acre  in  1984, 
and  half  of  these  are  greater  than  16  inches  d.b.h.,  only  50  trees  will  actually  be 
removed,  leaving  a  residual  stand  with  150  trees  per  acre. 

In  version  5.0,  the  SHRUB  and  COVER  extensions  have  been  functionally 
combined,  and  access  to  either  is  provided  by  inserting  the  COVER  keyword  in 
the  runstream  for  a  projection.  Keywords  associated  with  COVER  are 
described  by  Moeur  (1985). 

COVER  field  1:    Cycle  to  begin  COVER  calculations.  Once  begun, 

COVER  calculations  will  be  made  each  cycle  until  the 
end  of  the  projection;  default  =  beginning  of  projection, 
field  2:    Dataset  reference  number  for  COVER  output  file; 
default  =  18. 
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In  version  4.0,  probability  of  occurrence,  cover,  and  height  predictions  are 
made  for  16  shrub  species  on  three  habitat  types  (grand  iirlClintonia,  western 
redcedar'C/iVj^o/j/'a,  and  western  hemlocklClintonia).  In  version  5.0,  these  equa- 
tions have  been  replaced  and  expanded  to  represent  31  shrub  species  on  34  hab- 
itat types.  The  conifer  crown  width  and  foliage  models  are  identical  in  versions 
4.0  and  5.0. 

Interpreting  Stand  There  have  been  two  minor  changes  in  the  interpretation  of  data  on  the 

Information  Data         STDINFO  record.  In  version  4.0,  when  the  habitat  type  code  indicates  a  phase 

that  is  not  recognized  by  the  model,  the  habitat  code  is  ignored  (model  coeffi- 
cients for  type  999  are  used).  In  version  5.0,  phase  designations  are  ignored; 
the  code  is  rounded  to  habitat  type  and  coefficients  for  the  type  are  used. 

The  second  change  involves  assignment  of  default  values  when  more  than  one 
STDINFO  record  is  found  in  the  set  of  keyword  records  for  a  single  projection. 
In  version  4.0,  the  first  record  is  entirely  ignored  if  a  second  is  encountered. 
Thus,  blank  fields  are  assigned  the  default  parameter  values  indicated  in  the 
User's  Guide.  In  version  5.0,  blank  fields  on  second  and  subsequent  records  are 
assigned  values  from  data  entered  on  the  first  record. 
The  field  definitions  for  the  STDINFO  record  have  not  changed. 

Parameters  have  been  added  to  fields  4  and  5  of  the  MORTMULT  record  in 
order  to  specify  a  range  of  d.b.h.'s  to  which  the  multipliers  will  be  applied.  As 
in  version  4.0,  once  implemented,  the  multipliers  and  diameter  limits  are  in 
effect  for  remaining  cycles  in  a  projection  or  until  overridden  by  a  subsequent 
MORTMULT  record.  Also,  only  one  multiplier  and  one  set  of  diameter  limits 
can  be  specified  for  each  species  in  any  cycle. 

MORTMULT  field  1:    Cycle  in  which  mortality  multiplier  is  to  be  applied. 

Once  multipliers  take  effect,  they  remain  in  effect  until 
replaced  with  a  subsequent  request.  If  blank,  multipliers 
take  effect  at  the  start  of  the  projection. 

field  2:    Species  number  (see  appendix  A,  table  14)  to  which  mul- 
tiplier is  to  be  applied;  default  =  all  species, 
field  3:    The  value  of  the  multiplier  to  be  used;  default  =  1.0. 
field  4:    Minimum  d.b.h.  to  which  multiplier  will  be  applied; 
default  =  0.0. 

field  5:    Maximum  d.b.h.  to  which  multiplier  will  be  applied; 
default  =  999.0. 

Tree  Records  In  version  5.0,  the  interpretation  of  increment  data  has  changed  slightly. 

With  calibration  options  1  and  3  (see  Wykoff  and  others  1982,  p.  20),  incre- 
ments are  input  as  the  difference  between  two  successive  d.b.h.  or  height  meas- 
urements; the  value  entered  as  d.b.h.  or  height  increment  is  actually  a  d.b.h.  or 
height— the  program  automatically  performs  the  subtraction.  In  version  4.0 
with  calibration  options  1  and  3,  input  increments  of  0  are  treated  as  real  meas- 
urements and  input  increments  that  are  exactly  equal  to  the  d.b.h.  or  height 
measurements  are  ignored.  In  version  5.0,  only  non-zero  entries  in  the  incre- 
ment fields  are  treated  as  growth  samples  regardless  of  the  method  used  to 
enter  increment  data.  Further,  when  increment  is  computed  as  the  difference 
between  two  successive  height  or  d.b.h.  measurements,  all  negative  values  are 
rejected. 


Mortality 
Multipliers 
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Point-specific  site  descriptors  have  been  added  to  the  tree  records  as  an 
optional  input  to  facihtate  the  regeneration  estabUshment  model.  Data  that 
may  be  entered  include  slope,  aspect,  habitat  type,  topographic  position,  and 
site  preparation;  the  coding  of  these  variables  is  described  by  Ferguson  and 
Crookston  (1984,  p.  21).  The  default  tree  record  format  has  not  changed  other 
than  to  accommodate  the  optional  point-specific  site  descriptors. 

The  optional  input  is  read  from  the  tree  records  when  field  2  of  the 
TREEDATA  record  is  not  blank. 

TREE  DATA  field  1:    Dataset  reference  number  for  tree  record  input  file; 

default  =  2. 

field  2:  Any  numeric  value  in  this  field  will  cause  point-specific 
site  descriptors  to  be  read  from  tree  records.  Default  = 
point-specific  site  descriptors  are  not  read. 

NEW  OPTIONS 

Options  have  been  included  in  version  5.0  that  add  flexibiUty  to  activity 
schedules,  provide  better  description  of  the  input  data,  permit  variable  mer- 
chantability standards  for  computing  board  foot  volume,  and  afford  greater 
control  of  growth  functions. 

Event  Monitor  The  event  monitor  schedules  management  activities  conditional  on  the  status 

of  selected  stand  attributes  or  the  occurrence  of  certain  management  activities. 
The  IF  record  indicates  that  subsequent  records  contain  information  for  the 
event  monitor.  IF  is  followed  by  a  logical  expression  coded  on  one  or  more  sup- 
plemental data  records.  A  THEN  record  signals  the  end  of  the  logical  expres- 
sion and  is  followed  by  one  or  more  keywords  for  Prognosis  Model  activities. 
An  ENDIF  record  signals  the  end  of  event  monitor  information.  For  example, 

I  F 

(BBA      GT  120) 
THEN 

THINBBA  16  80 

ESTAB 

END 

END  I  F 

will  result  in  a  thinning  from  below  to  80  ft^  of  basal  area  followed  by  a  call  to 
the  regeneration  establishment  model  each  time  stand  basal  area  reaches 
120  ft2. 

Any  activity  that  can  be  scheduled  normally  can  be  alternatively  scheduled 
with  the  event  monitor;  the  scheduling  date  (field  1)  on  the  activity  keyword  is 
then  interpreted  as  a  delay.  If  the  date  entered  is  10,  the  activity  will  be  sched- 
uled 10  years  following  the  first  time  that  the  logical  expression  is  true. 

Crookston  (1985)  gives  further  examples  along  with  complete  instructions  for 
coding  keywords  and  logical  expressions. 

The  BFVOLUME  record  changes  the  merchantability  standards  for  the  board 
foot  volume  equations.  The  parameters  for  BFVOLUME  are  the  same  as  for 
VOLUME  in  version  4.0,  except  that  values  entered  for  stump  height  are 
ignored.  Minimum  d.b.h.  and  minimum  top  diameter  may  be  changed  by  cycle 
for  any  or  all  species. 

BFVOLUME  field  1:    Number  of  the  cycle  in  which  new  merchantabihty  stan- 
dards are  to  take  effect;  default  is  beginning  of 
projection. 


Variable 
Merchantability 
Standards  for  the 
Board  Foot  Volume 
Equations 
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field 


field 


field 


4 


3 


2 


Species  number  (see  appendix  A.  table  14)  for  the  spe- 
cies that  is  to  be  affected  by  the  merchantability  limits; 
default  is  all  species. 

Minimum  merchantable  d.b.h.  (inches).  Trees  with 
smaller  d.b.h.  will  have  no  board  foot  volume.  If  the 
number  entered  here  is  less  than  the  top  diameter  (field 
4),  the  value  specified  for  minimum  top  diameter  will  be 
used  for  minimum  d.b.h.  as  well;  default  =  6.0  for 
lodgepole  pine,  7.0  for  £J1  other  species. 
The  minimum  top  diameter  (inches);  default  =  4.5. 


Descriptive 
Statistics  for  Input 
Data 


Starting  Projections 
from  Bare  Ground 


Compressing  the 
Tree  List 


An  optional  output  table  gives  a  statistical  description  of  the  input  data  for  a 
projection.  This  table  is  requested  with  the  STATS  record.  Statistics  reported 
include  board  foot  volume,  cubic  foot  volume,  trees  per  acre,  and  basal  area  per 
acre  for  each  species  present  in  the  stand.  Also  given  are  the  mean,  standard 
deviation,  coefficient  of  variation,  and  confidence  Umits  across  sample  plots  for 
stand  totals  of  these  volume  and  density  measures.  Significance  level  is  input 
in  field  1  of  the  STATS  record;  the  corresponding  value  of  Student's  t  is  com- 
puted by  the  program. 

STATS  field  1:    Significance  level  for  computing  confidence  limits; 


With  the  Establishment  Model  in  place,  it  is  now  possible  to  begin  projec- 
tions without  tree  records.  To  do  this,  a  NOTREES  record  must  be  inserted  in 
the  runstream.  There  are  no  associated  parameters.  The  establishment  model 
(Ferguson  and  Crookston  1984)  can  then  be  used  to  generate  a  list  of  seedlings 
based  on  stand  parameters  and  proposed  treatment  strategies. 

When  there  are  many  records  in  the  tree  list,  the  record-compression 
algorithm  can  reduce  the  number  of  records  while  minimizing  loss  of  within- 
stand  variation  in  tree  attributes  that  are  important  for  increment  prediction. 
As  a  result,  the  number  of  calculations  required  for  the  projection  (and  the  cost 
of  the  projection)  can  be  substantially  reduced  without  significant  bias.  Com- 
pression of  the  tree  list  is  requested  with  the  COMPRESS  record;  there  are 
four  associated  parameters: 

COMPRESS    field  1:    Cycle  in  which  tree  list  will  be  compressed;  default  is 

the  beginning  of  the  projection. 

field  2:    Number  of  tree  records  that  will  remain  following  com- 
pression; default  =  150.0. 

field  3:    Percentage  of  new  records  that  will  be  determined  by 
finding  the  largest  gaps  in  the  classification  space. 
Remainder  of  records  will  result  from  splitting  the 
classes  with  the  greatest  variation;  default  =  50.0. 

field  4:    Any  numeric  entry  will  cause  debug  output  to  be 

printed  for  the  compression  algorithm;  default  is  no 
compression  debug  output. 

If  COMPRESS  is  specified  for  the  first  or  second  cycle,  the  tree  records  will 
be  subsequently  tripled  unless  NOTRIPLE  or  NUMTRIP  (see  Wykoff  and 
others  1982,  p.  93)  is  used  to  suppress  the  record  tripling  feature. 
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Simulating  Top  Two  options  were  included  in  version  5.0  to  simulate  the  effects  of  defoliation 

Mortality  by  insects  when  interactive  population  dynamics  models  are  unavailable.  By 

using  the  HTGSTOP  record,  height  growth  may  be  set  to  zero  for  randomly 
selected  tree  records  with  heights  between  user-specified  limits.  The  TOPKILL 
record  provides  for  permanently  removing  a  randomly  selected  proportion  of 
the  top  as  well.  There  are  five  parameters  on  the  HTGSTOP  record  to  indicate 
timing,  species  affected,  upper  and  lower  height  limits,  and  the  probabiUty  that 
the  tree  will  not  grow  in  height.  These  same  five  parameters  are  entered  on  the 
TOPKILL  record  along  with  a  mean  and  standard  deviation  for  the  distribution 
describing  the  proportion  of  total  tree  height  lost  to  TOPKILL.  This  distribu- 
tion is  assumed  to  be  Normal,  but  in  no  case  will  a  tree's  crown  ratio  be 
reduced  to  less  than  5  percent  or  wiU  the  maximum  change  in  total  height  be 
greater  than  25  percent. 

HTGSTOP 

TOPKILL       field  1:    Cycle  in  which  top-kill  and/or  growth  loss  wiU  be  simu- 
lated; default  is  first  cycle. 

field  2:    Species  that  will  be  affected;  default  is  all  species. 

field  3:    Lower  Umit  of  range  of  heights  that  wiU  be  affected; 
default  =  0.0. 

field  4:    Upper  Hmit  of  range  of  heights  that  wiU  be  affected; 
default  =  0.0. 

field  5:    Probability  that  a  tree  will  experience  top-kill  or  growth 

loss;  default  =  0.0. 
For  TOPKILL  record  only: 

field  6:    Mean  of  the  distribution  of  the  proportion  of  total  tree 

height  lost  to  top-kill;  default  =  0.0. 
field  7:    Standard  deviation  of  the  distribution  of  the  proportion 

of  total  tree  height  lost  to  top-kill;  default  =  0.0. 

Top  damage  will  only  be  simulated  in  the  cycle  that  is  specified  on  the 
HTGSTOP  or  TOPKILL  record.  Periodic  infestations  can  be  simulated  by 
including  HTGSTOP  or  TOPKILL  records  for  additional  cycles. 

CHANGES  IN  PROGRAM  OUTPUT 

Although  the  format  of  the  base  model  output  has  not  changed  greatly  in 
version  5.0,  additions  to  the  base  model  have  resulted  in  new  tables.  There 
have  also  been  subtle  changes  to  the  output  that  result  from  new  or  expanded 
options  £md  from  changes  in  the  interpretation  of  input  data  on  the  frequency 
of  nonstockable  plots. 

The  regeneration  estabhshment  model  and  the  SHRUBS  and  CANOPY 
models  produce  output  tables  that  summarize  the  operation  of  these  models 
when  they  are  invoked.  These  tables  are  explained  in  the  appropriate  User's 
Manual  (Ferguson  and  Crookston  1984;  Moeur  1985).  The  STATS  option 
produces  another  new  table  that  displays  the  distribution  of  stand  attributes 
by  species  and  by  sample  plot  when  more  than  one  plot  is  used  to  inventory 
the  stand  (fig.  5). 

The  change  in  the  interpretation  of  input  data  on  the  frequency  of  non- 
stockable plots  resulted  in  adjustments  to  the  per-acre  values  reported  for 
st£md  attributes.  As  was  earlier  described,  trees  per  acre  represented  by  a  tree 
record  is  now  computed  on  the  basis  of  stockable  area.  When  nonstockable 
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SPEC  I £S 


GENLRAL  SPECIES  SUMMARY   FOR   THE  CRUISE 
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Figure  5.  — Prognosis  Model  output  table  produced  by  the  STATS  option. 

points  are  present  in  an  inventory,  the  "STAND  COMPOSITION"  table 
(Wykoff  and  others  1982;  fig.  7.  pp.  34-35)  and  the  table  displaying  "ATTRIB- 
UTES OF  SELECTED  SAMPLE  TREES"  and  "ADDITIONAL  STAND 
ATTRIBUTES"  (Wykoff  and  others  1982;  fig.  8,  pp.  36-37)  now  display  per- 
acre  stand  and  tree  attributes  on  the  basis  of  stockable  area.  The  table  head- 
ings have  been  modified  to  reflect  the  change  (figs.  6  and  7).  The  "SUMMARY 
STATISTICS"  table  (Wykoff  and  others  1982;  fig.  9.  p.  37)  still  reports  per-acre 
stand  attributes  on  the  basis  of  toted  stand  area.  The  table  heading,  however, 
has  been  modified  for  contrast  (fig.  8).  The  trees  per  acre  reported  for  each  tree 
record  when  the  TREELIST  option  is  used  is  still  based  on  total  stand  area  as 
well. 

If  nonstockable  plots  are  present  in  an  inventory,  the  "OPTIONS 
SELECTED  BY  DEFAULT"  table  now  reports  the  ratio  of  stockable  to  total 
plots  and  the  basis  for  computing  per-acre  attributes  in  each  of  the  major 
tables  (fig.  9). 

The  final  change  in  program  output  results  from  the  i;nplementation  of  the 
regeneration  establishment  model,  the  compression  aJgorithm,  and  the  complete 
removal  of  tree  records  in  thinning.  When  tree  records  are  added  or  deleted,  a 
realignment  of  the  records  displayed  in  the  "ATTRIBUTES  OF  SELECTED 
SAMPLE  TREES  '  table  is  necessary.  Changes  are  noted  by  placing  asterisks 
(**)  following  the  date  that  records  are  reselected.  An  explematory  note  is  also 
appended  to  the  table  (fig.  7). 

CONCLUSIONS 

There  are  si'.bstantial  changes  in  program  operation  and  model  formulation 
represented  in  version  5.0.  As  a  result,  the  program  affords  both  greater  flexi- 
bility and  improved  performance.  Permanent  sample  plots  were  used  to  com- 
pare version  5.0  with  previously  available  methods  of  growth  projection  for 
Northern  Rocky  Mounteiin  conifers.  Results  show  that  version  5.0  projections 
of  stand  attributes  have  smaller  biases  and  smaller  error  variemces  thgin  do 
projections  produced  by  other  methods  (Wykoff  1985;  Stage  in  preparation). 
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STAND   ID:  S248112 


STAND  GROWTH   PROGNOSIS  SYSTEM  VERSION  5.1  INLAND  EMPIRE 

MANAGEMENT  CODE:    NONE  SHELTERWOOD   PRESCRIPTION   FROM  THE  USER'S  MANUAL 


STAND  COMPOSITION   (BASED  ON   STOCKABLE  AREA) 
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DISTRIBUTION  OF  STAND  ATTRIBUTES  BY  DBH 


10 


30 


50  70 
(DBH    IN  INCHES) 


90 


100 


TOTAL/ACRE 

OF  STAND 
ATTR I BUTES 


DISTRIBUTION  OF  STAND  ATTRIBUTES  BY 
SPECIES  AND  3   USER-DEFINED  SUBCLASSES 


1977  TREES 

0. 

.  1 

0, 

,  1 

3, 

.2 

6, 

.  1 

8 

.5 

12. 

.7 

590. 

TREES 

VOLUME: 

TOTAL 

5. 

,8 

7. 

.9 

9. 

10. 

,0 

1 1 

.5 

12. 

,7 

1623. 

CUFT 

MERCH 

8. 

2 

9. 

9. 

6 

10. 

.9 

1  1 

.5 

12. 

7 

1118. 

CUFT 

MERCH 

8. 

0 

9. 

4 

9. 

,6 

10. 

9 

1 1 

.5 

12. 

7 

l*57U. 

BDFT 

REMOVAL 

0 

.  1 

0. 

,  1 

0. 

,  1 

1 , 

.  2 

3, 

.2 

10, 

.4 

326. 

TREES 

VOLUME: 

TOTAL 

8, 

.0 

9. 

6 

9. 

6 

10, 

.1* 

10, 

.1* 

10. 

,1* 

255. 

CUFT 

MERCH 

9, 

.6 

9. 

6 

9. 

6 

10. 

.1* 

10, 

.1* 

10. 

,1* 

217. 

CUFT 

MERCH 

8. 

.0 

9. 

6 

9. 

6 

10. 

,4 

10. 

,4 

10. 

,4 

911  . 

BDFT 

RES  I  DUAL 


ACCRETION 
MORTALITY 


4.0 


5.3 
6.  1 


5.3 


6.2 
7.9 


6.2 


7.9 
8.5 


7.9 


9.5 


9.4  10.9 

9.5  11.5 


12.7 


12.7 
12.7 


264.  TREES 


27 

■  % 

DF2, 

15 

.% 

GF2, 

15, 

■  % 

GFl  , 

10 

24, 

■  % 

LP1  . 

21 

.% 

GFl  , 

20, 

.% 

DFl  , 

12 

32. 

■  % 

LPl  , 

24, 

.% 

DF1  , 

15, 

,% 

LI  , 

10 

31  . 

■  % 

LPl  , 

23, 

.% 

DFl  , 

17, 

.% 

LI  , 

10 

48. 

.% 

DF2, 

28. 

.% 

GF2, 

18. 

■  % 

C2, 

4 

47. 

% 

LP2, 

40. 

% 

DF2, 

9. 

% 

L2, 

4 

51  . 

% 

LP2, 

41  . 

% 

DF2, 

8. 

% 

L2, 

0 

50. 

% 

LP2, 

37. 

% 

DF2, 

13. 

% 

L2, 

0 

33. 

% 

GF1  , 

19. 

% 

CI  , 

19. 

% 

DFl  , 

13 

31  . 

% 

GFl  , 

25. 

% 

DFl  , 

13. 

% 

LPl  , 

1  1 

39. 

% 

LPl  , 

25. 

% 

LI  . 

20. 

% 

GFl  , 

1  3 

%  C2 

%  LI 

%  LP2 

%  LP2 

%  L2 


1987  TREES 

5, 

.6 

6, 

.8 

7, 

.9 

9 

.4 

1 1 

.  2 

16, 

,  2 

VOLUME: 

TOTAL 

7. 

,  3 

8, 

,9 

9, 

.9 

1 1  , 

.6 

13 

.7 

16, 

,  2 

MERCH 

7. 

,  7 

9. 

,4 

10. 

,4 

12, 

.0 

13, 

.  7 

16. 

,  2 

MERCH 

7. 

.9 

9. 

,4 

10. 

,6 

12, 

.4 

13, 

.9 

16. 

,2 

ACCRETION 

6. 

6 

7. 

7 

9. 

0 

1  1  . 

,0 

13. 

,  1 

16. 

2 

MORTALITY 

6. 

8 

8. 

9 

9. 

4 

10. 

6 

12. 

,5 

16. 

2 

242. 

TREES 

33. 

.% 

GFl  , 

20, 

.% 

CI  , 

19. 

% 

DFl  , 

12. 

,% 

LI 

2145. 

CUFT 

28. 

% 

GFl  , 

24. 

.% 

DFl  , 

21  . 

% 

LPl  , 

12. 

■  % 

LI 

1763. 

CUFT 

26. 

% 

GFl  , 

25. 

,% 

DFl  , 

24. 

% 

LPl  , 

13. 

■  % 

LI 

7247. 

BDFT 

26. 

% 

LPl  , 

25. 

■  % 

DFl , 

24. 

% 

GFl  , 

15. 

% 

LI 

116. 

CUFT/YR 

36. 

% 

GFl , 

22. 

% 

DFl  , 

12. 

% 

CI  , 

12. 

% 

WH1 

18. 

CUFT/YR 

36. 

% 

LPl , 

23. 

% 

GFl  , 

21  . 

% 

LI  , 

15. 

% 

DFl 

1997  TREES 

7, 

,2 

8, 

,4 

9, 

.6 

1 1  , 

.  1 

13. 

,4 

19, 

,2 

222. 

TREES 

33. 

% 

GFl  , 

21  . 

% 

CI  , 

19. 

.% 

DFl , 

1  1  , 

.% 

VOLUME: 

TOTAL 

8, 

,4 

9. 

,9 

1 1  , 

.2 

13, 

,4 

15. 

5 

19. 

,2 

3123. 

CUFT 

31  . 

% 

GFl  , 

24. 

% 

DFl  , 

16. 

% 

LPl  , 

10, 

,% 

MERCH 

8, 

,6 

10. 

,  2 

1 1  , 

.5 

13, 

.4 

15. 

6 

19. 

,2 

2826. 

CUFT 

31  . 

% 

GFl  , 

23. 

% 

DFl  , 

17. 

■  % 

LPl  , 

1  1  , 

,% 

MERCH 

9. 

,  2 

10. 

,5 

1 1 . 

,  7 

13, 

.5 

15. 

8 

19. 

,  2 

1 1998. 

BDFT 

30. 

% 

GFl , 

25. 

% 

DFl  , 

18. 

1 

LPl  , 

1  1  , 

■  % 

ACCRETION 

8. 

,0 

9. 

6 

1 1  . 

,  1 

12. 

,5 

15. 

1 

19. 

2 

141  . 

CUFT/YR 

41  . 

% 

GFl  , 

20. 

% 

DFl  , 

14. 

% 

WH1  , 

13. 

■  % 

MORTALITY 

8. 

,0 

9. 

8 

10. 

8 

12. 

,  0 

14. 

3 

19. 

2 

28. 

GUFT/YR 

30. 

% 

LPl  , 

24. 

% 

GFl  , 

20. 

% 

LI  , 

18. 

■  % 

2007 


TREES 
VOLUME: 
TOTAL 
MERCH 
MERCH 


8. 

,2 

9 

.9 

1 1  . 

,  3 

12, 

,  5 

15, 

.4 

21 

.  7 

201  . 

TREES 

34. 

%  GFl , 

23. 

%     CI  , 

19. 

.% 

DFl  . 

10, 

.% 

9. 

,6 

1  1  , 

.6 

12. 

,8 

15. 

,0 

17, 

.  1 

21 

.  7 

4248. 

CUFT 

35. 

%  GFl, 

23. 

%  DFl, 

12. 

% 

LPl  , 

1  1  , 

.% 

9. 

,9 

1 1  , 

.  7 

13. 

3 

15. 

,  1 

17, 

,  2 

21 

.7 

3987. 

CUFT 

35. 

%  GFl, 

23. 

%  DFl, 

12. 

% 

LPl  , 

1  1  , 

■  % 

10. 

.0 

1 1  , 

.8 

13. 

4 

15. 

,4 

17. 

,4 

21  , 

.7 

17703. 

BDFT 

34. 

%  GFl, 

24. 

%  DFl  , 

12. 

% 

LPl  , 

1  1  . 

,% 

LI 

CI 


REMOVAL 
VOLUME: 
TOTAL 
MERCH 
MERCH 


ACCRET ION 
MORTALITY 


7, 

.2 

10, 

.0 

1 1 

.2 

1 1  . 

,8 

13, 

.4 

21 

.6 

44. 

TREES 

43. 

.% 

LI  , 

32, 

.% 

LPl  , 

1  1  . 

% 

CI  , 

10, 

,% 

10. 

,  0 

1 1  , 

.2 

1 1 

.8 

12. 

8 

14. 

.5 

21 

.6 

897. 

CUFT 

55. 

■  % 

LPl  , 

39, 

■  % 

LI  , 

3. 

% 

DFl  , 

2. 

.% 

10. 

.  1 

1  1  , 

.  3 

1 1 

.8 

12. 

8 

14, 

.5 

21 

.6 

828. 

CUFT 

57. 

■  % 

LPl  , 

40, 

■  % 

LI  , 

2. 

% 

CI  , 

1  . 

■  % 

10. 

,  3 

1  1  , 

.  3 

1 1  , 

.8 

13. 

1 

14. 

,  7 

21 

.6 

3673. 

BDFT 

57. 

■  % 

LPl  , 

40, 

■  % 

LI  , 

1  . 

% 

CI  , 

1  . 

■  % 

8. 

2 

9. 

,  7 

1  1  , 

.  3 

13. 

3 

15. 

,  7 

21 

.  7 

157. 

TREES 

42. 

% 

GFl  , 

26. 

■  % 

CI  , 

22. 

% 

DFl  , 

10. 

% 

9. 

,4 

1  1  . 

,  3 

12, 

.6 

14. 

5 

17. 

0 

21  , 

.  7 

144. 

CUFT/YR 

42. 

% 

GFl  , 

20. 

■  % 

WH1  , 

19. 

% 

CI  , 

18. 

% 

9. 

,0 

1 1  . 

,  3 

12. 

,  7 

15. 

1 

17. 

,2 

21  , 

,  7 

19. 

CUFT/YR 

46. 

% 

GFl  , 

34, 

■  % 

DFl , 

9. 

% 

CI  , 

9. 

% 

2017  TREES 

9. 

,6 

1 1  , 

,4 

13. 

,  1 

15, 

,  1 

17, 

.  8 

23. 

,  3 

147. 

TREES 

42. 

.% 

GFl , 

26. 

.% 

CI  , 

21 

.% 

DFl  , 

10. 

■  % 

WHl 

VOLUME: 

TOTAL 

1  1  . 

,  1 

13. 

,  1 

15. 

,  7 

17, 

,  3 

19, 

.2 

23. 

3 

4606. 

CUFT 

43. 

.% 

GFl  , 

25. 

■  % 

DFl  , 

16, 

.% 

WHl  , 

16. 

% 

CI 

MERCH 

1  1  . 

1 

13. 

,  3 

15. 

9 

17. 

,  3 

19, 

,  2 

23. 

3 

4397. 

CUFT 

43. 

.% 

GFl  , 

25. 

■  % 

DFl  , 

16, 

■  % 

WH1  , 

15. 

% 

CI 

MERCH 

1  1  . 

4 

13. 

,7 

16. 

3 

17. 

,6 

20. 

,0 

23. 

3 

20929. 

BDFT 

43. 

■  % 

GFl  , 

26. 

■  % 

DFl  , 

17, 

.% 

WHl  , 

13. 

% 

CI 

ACCRET ION 

10. 

8 

12. 

,9 

14. 

8 

16. 

.5 

19. 

,2 

23. 

3 

169. 

CUFT/YR 

44. 

■  % 

GFl  , 

23. 

% 

WH1  , 

17, 

.% 

CI  . 

15. 

% 

DFl 

MORTALITY 

9. 

9 

12. 

8 

14. 

4 

16. 

5 

18. 

,  7 

23. 

3 

32. 

CUFT/YR 

50. 

■  % 

GFl  , 

30. 

% 

DFl  , 

9, 

■  % 

WHl  , 

9. 

% 

CI 

Figure  6.— Stand  composition  table  from  the  Prognosis  Model  output.  This  table 
was  produced  by  the  same  set  of  keyword  and  tree  records  that  were  used  for  an 
example  by  Wykoff  and  others  (1982;  fig.  7,  pp.  34-35).  Note  the  change  in  table 
heading  that  reflects  the  change  in  treatment  of  stockability  data.  Note  also  that  a 
nonstockable  plot  in  the  inventory  resulted  in  a  change  in  initial  conditions. 
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Figure  7.— Tree  and  stand  attributes  table  fronn  the  Prognosis  Model  output.  As  for 
figure  6,  input  conditions  nnatcti  those  used  to  produce  tables  in  the  example  by 
Wyl<off  and  others  (1982;  fig.  8,  pp.  36-37).  Note  that  the  heading  has  changed  to 
Show  that  per-acre  values  are  based  on  stockable  area  and  that  initial  conditions 
are  somewhat  different.  Note  also  that  when  the  regeneration  establishment  model, 
the  compression  algorithm,  and/or  thinning  result  in  the  addition  or  removal  of  tree 
records,  new  tree  records  are  selected  for  display.  The  year  is  flagged  with  two 
asterisks  (see  the  year  2047)  and  an  explanatory  note  is  appended  to  the  table. 
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Figure  7  —  (con.) 
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Figure  8  — Summary  statistics  table  from  the  Prognosis  Model  output.  Input  condi- 
tions are  the  same  as  displayed  in  Wykoff  and  others  (1982:  fig.  9:  p.  37).  The  table 
heading  was  modified  to  state  that  statistics  are  based  on  total  stand  area. 


OPTIONS  SFi rcTEn  BY  nrFAUl  r 


TREEFMT         (23X,H4,3X,    F2.0,I1,    A3  ,  F  3  .  1  ,  F  2  .  1  ,  3X  ,  F  3  .  0  .  T6  3  ,  F  3  .  0       .  r6() ,  F  3  .  1  ,  1 1|8  .      II, 3X.  12, 
2  11,  T66, 2  11,13,  211) 

DESIGN  BASAI    AREA   FACTOR^   UO.O;    INVERSE   OF    FIXED   PLOT   AREA--      300.0;    BREAK  DBH-^  5.0 

NUMBER  OF    PLOTS=  11;    NGN- S I OCKABL E    PLOTS-  1;    STAND  SAMPLING  WEIGHT^  11.00000 

STAND  ATTRIBUTES  ARE  CAIXUIAIED   PER  ACRE  Of    STOCKABLE  AREA.      STAND  STATISTICS 
IN  SUMMARY  TABLE  ARE  MULTIPLIED  BY  0.909   TO    INCLUDE   TOTAL  STAND  AREA. 


Figure  9.  — Prognosis  Model  output  table  showing  options  selected  by  default.  This 
figure  illustrates  the  message  produced  when  nonstockable  plots  are  encountered 
in  the  inventory.  Contrast  to  Wykoff  and  others  (1982;  fig.  6,  pp.  32-33). 
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APPENDIX  A:  SUMMARY  OF  CODES  USED  IN  THE 
PROGNOSIS  MODEL 


Table  12.— Codes  for  the  Forests  represented  in  the  Inland 
Empire  version  of  the  Prognosis  Model 


Forest 

Code 

Forest 

Code 

Bitterroot 

3 

Kanlksu 

13 

Idaho  Panhandle^ 

4 

Kootenai 

14 

Clearwater 

5 

Lolo 

16 

Coeur  d'Alene 

6 

Nezperce 

17 

Colville 

7 

St.  Joe 

18 

Flathead 

10 

'The  Idaho  Panhandle  National  Forests  is  an  administrative  unit 
that  encompasses  the  Coeur  d'Alene.  Kaniksu,  and  St.  Joe  National 
Forests.  All  of  the  growth  models  use  coefficients  for  the  St.  Joe 
National  Forest  when  the  Panhandle  is  specified. 


Table  13.— Codes  for  habitat  types  represented  in  the  Inland  Empire 
version  of  the  Prognosis  Model^ 


Code 

Abbreviation 

Habitat  type  name 

130 

PIPO/AGSP 

Pinus  ponderosa/Agropyron  spicatum 

170 

PIPO/SYAL 

Pinus  ponderosa/Symphohcarpos  albus 

250 

PSfvlE/VACA 

Pseudotsuga  menziesii/Vaccinium  caespitosum 

260 

PSME/PHMA 

Pseudotsuga  menziesii/Physocarpus  malvaceus 

280 

PSME/VAGL 

Pseudotsuga  menziesii/Vaccinium  globulare 

290 

PSME/LIBO 

Pseudotsuga  menziesii/Linnaea  borealis 

310 

PSME/SYAL 

Pseudotsuga  menziesii/Symphoricarpos  albus 

320 

PSME/CARU 

Pseudotsuga  menziesli/Calamagrostis  rubescens 

330 

PSIVIE/CAGE 

Pseudotsuga  menziesii/Carex  geyeri 

420 

PICEA/CLUN 

Picea/Clintonia  uniflora 

470 

PICEA/LIBO 

Picea/Linnaea  borealis 

510 

ABGR/XETE 

Abies  g rand is/Xeroplny Hum  tenax 

520 

ABGR/CLUN 

Abies  grandis/Clintonia  uniflora 

530 

THPL/CLUN 

Thuja  plicata/Clintonia  uniflora 

540 

THPL/ATFI 

Thuja  plicata/Athyrium  filix-femina 

550 

THPL/OPHO 

Thuja  plicata/Oplopanax  horridum 

570 

TSHE/CLUN 

Tsuga  heterophylla/Clintonia  uniflora 

610 

ABLA/OPHO 

Abies  lasiocarpa/Oplopanax  horridum 

620 

ABLA/CLUN 

Abies  lasiocarpa/Clintonia  uniflora 

640 

ABLA/VACA 

Abies  lasiocarpa/Vaccinium  caespitosum 

660 

ABLA/LIBO 

Abies  lasiocarpa/Linnaea  borealis 

670 

ABLA/MEFE 

Abies  lasiocarpa/Menziesia  ferruginea 

680 

TSIVIE/MEFE 

Tsuga  mertensiana/t^enziesia  ferruginea 

690 

ABLA/XETE 

Abies  lasiocarpa/Xerophyllum  tenax 

710 

TSf^E/XETE 

Tsuga  mertensiana/Xerophyllum  tenax 

720 

ABLA/VAGL 

Abies  lasiocarpa/Vaccinium  globulare 

730 

ABLA/VASC 

Abies  lasiocarpa/Vaccinium  scoparium 

830 

ABLA/LUHI 

Abies  lasiocarpa/Luzula  hitchcockii 

850 

PIALABLA 

Pinus  albicaulis/Abies  lasiocarpa 

999 

OTHER 

'From  Pfister  and  others  (1977). 
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Table  14.— Tree  species  recognized  by  the  Prognosis  Model  with 
default  coding  conventions 


Default 

Numeric 

Common  name 

Scientific  name 

input  code 

code 

Western  white  pine 

Pinus  monticola 

WP 

1 

Western  larch 

Larix  occidentalis 

L 

2 

Douglas-fir 

Pseudotsuga  menziesii 

DF 

3 

Grand  fir 

Abies  grandis 

GF 

4 

Western  hemlock 

Tsuga  heterophylla 

WH 

5 

Western  redcedar 

Thuja  pi  i  cat  a 

c 

6 

Lodgepole  pine 

Pinus  contorta 

LP 

7 

Engelmann  spruce 

Picea  engelmannii 

S 

8 

Subalpine  fir 

Abies  lasiocarpa 

AF 

9 

Ponderosa  pine 

Pinus  ponderosa 

PP 

10 

fvlountain  hemlock 

Tsuga  mertensiana 

11 

Table  15.— Aspect  codes  '•'s'^'®  16.— Slope  codes 


Aspect  Azimuth  (degrees)      Code  Slope  angle  (%)  Code 


North  337.5  -  22.5  1  <  -  5  0 
Northeast  22.6  -  67.5  2  6  -  15  1 
East  67.6  -  112.5  3  16  -  25  2 
Southeast  112.6  -  157.5  4  26  -  35  3 
South  157.6  -  202.5  5  36  -  45  4 
Southwest  202.6  -  247.5  6  46  -  55  5 
West  247.6  -  292.5  7  56  -  65  6 
Northwest  292.6  -  337.5  8  66  -  75  7 
Level  9  76  -  85  8 
  >     86  9 


Table  17.— Crown  ratio  codes 


Crown  ratio  (%) 

Code 

1  -  10 

1 

11-20 

2 

21  -  30 

3 

31  -  40 

4 

41  -  50 

5 

51  -  60 

6 

61  -  70 

7 

71  -  80 

8 

>  81 

9 

Table  18.- 

-Interpreting  damage  codes 

(IDCO) 

Code 

Interpretation 

73 

Tree  top  is  missing 

74 

Tree  top  is  dead 

all  others 

Ignored 

Table  19.— Interpreting  tree  history  codes  (ITH) 
Code  Interpretation 

5  Tree  died  during  mortality  observation  period;  record  is  used  to 

backdate  density  for  model  scaling,  but  is  not  projected. 

6,7  Tree  died  prior  to  mortality  observation  period;  record  is  ignored. 

9  Special  record  (planar  intercept  in  Northern  Region  inventory); 

record  is  ignored. 

1,2,3,4,8    Various  categories  of  live  trees;  records  are  projected. 
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Table  20.— Interpreting  tree  value  codes  (IMC) 


Code  Interpretation 


1  Desirable  tree 

2  Acceptable  tree 

3  Live  cull 

8  Nonstockable  point 


All  other  codes  are  interpreted  as  3 


APPENDIX  B:  SUMMARY  OF  KEYWORDS,  ASSOCIATED 
PARAMETERS,  AND  DEFAULT  CONDITIONS 

Note:  Appendix  B  contains  summaries  of  keywords  that  are  presented  in 
both  the  version  5.0  supplement  and  the  original  User's  Guide  (Wykoff  and 
others  1982).  For  a  more  detailed  description,  refer  to  the  page  number  given 
beneath  each  keyword.  Numbers  preceded  by  the  letter  "S"  are  for  either  new 
or  updated  options  and  the  indicated  pages  are  found  in  this  supplement. 
Otherwise,  refer  to  the  User's  Guide. 

Rules  for  Coding  1.  All  option  keywords  start  in  column  1. 

Keyword  Records  2.  The  numerical  values  (parameters)  needed  to  implement  an  option  are 

contained  in  seven  numeric  fields  that  are  10  columns  wide.  The  first  parameter 
field  begins  in  column  11.  Values  should  either  be  right-justified  in  the  numeric 
field  or  be  followed  by  a  decimed  point. 

3.  Blank  numeric  fields  are  not  treated  as  zeroes.  If  a  blank  field  is  found,  the 
default  value  will  be  used.  If  zeroes  are  to  be  specified,  they  must  be  punched. 
Thus,  only  the  numeric  values  that  are  different  from  the  default  parameter 
values  need  to  be  specified. 

4.  All  supplemental  data  records  associated  with  a  keyword  must  be  provided 
if  the  keyword  is  used. 

5.  When  two  or  more  conflicting  options  are  specified,  the  last  one  specified 
will  be  used. 


Controlling  Program  Execution 

Keyword 

(page  Default  parameter 

reference)  Keyword  use  and  associated  parameters  or  conditions 


INVYEAR  Specify  the  starting  date  for  a  projection. 

(8)  field  1:  Year  in  which  simulation  is  to  begin. 

NUMCYCLE  Specify  the  number  of  cycles  in  a  projection. 
(8)  field  1:  Number  of  cycles  to  be  projected; 

maximum  number  of  cycles  is  40. 
PROCESS  Marks  the  end  of  an  input  file  for  a  single  projection  in 

(8)  a  runstream  and  triggers  the  beginning  of  the  simulation. 

Must  be  present  or  projection  will  not  run. 

STOP  Signal  the  end  of  Prognosis  Model  runstream. 

(8) 

TIME  INT  Specify  the  length  of  any  or  all  projection  cycles. 

(8)  field  1:  Number  of  a  cycle  whose  length  is  to  be 

changed. 

field  2:  Number  of  years  to  be  simulated  in  the 
cycle(s)  referenced  in  field  1. 


Change  all  cycles 


10  years 
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Entering  Stand  and  Tree  Characteristics 

Keyword 

(page  Default  parameter 

reference)  Keyword  use  and  associated  parameters  or  conditions 


DESIGN 
(10,  S-15) 


GROWTH 
(20) 


MGMTID 
(11) 


SPCODES 
(19) 


STDIDENT 
(11) 


STDINFO 
(12,S-17) 


TREEDATA 
(18,  S-17) 


field 

2: 

field 

3: 

field 

4: 

field 

5: 

field 

6: 

Enter  inventory  design  parameters, 
field  1:  If  positive,  basal  area  factor  for  variable 
radius  plots;  if  negative,  inverse  of  fixed 
area  plot  for  sampling  large  trees. 
Inverse  of  fixed  plot  area  for  measuring 
small  trees. 

d.b.h.  separating  large  from  small  trees. 
Number  of  plots  used  to  inventory  stand. 
Number  of  nonstockable  plots  in  stand 
inventory. 

Stand  weight  for  aggregation  of  projections. 
Identify  methods  used  to  measure  and  input  mortality 
and  height  and  diameter  increment  data, 
field  1:  Method  used  to  measure  diameter  increment, 
field  2:  Length  of  diameter  increment  measurement 
period. 

field  3:  Method  used  to  measure  height  increment, 
field  4:  Length  of  height  increment  measurement 
period. 

field  5:  Length  of  mortality  observation  period. 
Enter  an  alphanumeric  code  to  identify  the  silvicultural 
treatment  simulated  in  a  projection.  The  code  does  not 
affect  the  projection  but  is  printed  with  each  output 
table  and  on  each  line  in  the  summary  table. 
Supplemental  record:  enter  management  identifier  in 
first  four  columns. 
Identify  species  codes  used  on  the  input  tree  records, 
field  1:  Numeric  code  for  the  species  for  which  the 

code  is  to  be  changed. 
Supplemental  record:  Species  code  or  codes,  left 
justified  in  consecutive  4-column  fields.  If  all  codes 
are  replaced,  they  must  be  entered  in  order  of  numeric 
code.  If  only  one  code  is  replaced,  it  is  entered  in 
the  first  4  columns. 
Enter  steind  identification  code  and  descriptive  title  to 
label  the  output. 
Supplemental  record:  Stand  identification  code  is 
entered  in  columns  1-8;  title  is  entered  in  columns  9-80. 
Enter  data  that  describe  the  site  on  which  stand  is 
located. 

field  1:  National  Forest  on  which  stand  is  located. 

field  2:  Stand  habitat  type  code. 

field  3:  Stand  age. 

field  4:  Stand  aspect  code. 

field  5:  Stand  slope  code. 

field  6:  Stand  elevation  code. 

field  7:  Stand  site  index. 
Read  tree  data  from  dataset  referenced  by  the  unit 
number  recorded  in  field  1. 

field  1:  Dataset  reference  number. 

field  2:  Indicate  that  point-specific  site  descriptors 
are  included  with  tree  records. 


40  ft^/acre/tree 

300  plots/acre 
5  inches 
Count  the  plots 

Count  the  plots 
Number  of  plots 


0  (past  increment) 
10  years 

0  (past  increment) 

5  years 
5  years 

Default  code  is  "NONE' 
(MGMTID  record  not 
input);  if  supplemental 
record  is  blank, 
management  identifiers 
not  printed 


Change  for  all  species. 
Default  vEilues  are  given 
in  table  14;  a  blank  entry 
on  the  supplemental  record 
will  be  interpreted  as  a 
blank 


18  (St.  Joe) 

260  (PSME/PHMA) 

0  years 

9  (level) 

0  (<5%) 

38  hundred  feet 

0 


No  point-specific 
site  descriptors 
on  tree  records 


(con. 
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Entering  Stand  and  Tree  Characteristics  (Con.) 


Keyword 

(page 
reference) 


Keyword  use  and  associated  parameters 


Default  parameter 
or  conditions 


TREEFMT  Provide  a  format  statement  that  describes  the  layout  of 
(19)  a  tree  record. 

Two  supplemental  records:  A  FORTRAN  execution 

time  format  statement. 


See  table  5,  Wykoff 
and  others  1982 


Modifying  Tree  Volume  Calculations 

Keyword 


(page 

Default  parameter 

reference) 

Keyword  use  and  associated  parameters 

or  conditions 

nr  r  ULilN 

Enter  species-specific  parameters  for  log-hnear  form  and 

defect  correction  e<juation  for  board  foot  volume  estimates 

(24) 

(BFFDLN)  or  merchantable  cubic  foot  volume  estimates 

(MLr  ULIN ). 

field  1:  Numeric  code  for  the  species  for  which  the 

cn.|UclLlUll  lo  tU  Uc  Clldil^cvi.    X  lie  UcldLUt  CljUclLlUll 

*?iinnlip<?  ri  mnltinlipr  of  1  0  for  p?if*)i  <?npcips 

ChanP'P  all  ^nprips 

\^  lltJUl^L^    CUA  i^L^^^A^O 

field  2:  Intercept  term  for  log'linear  eQuation. 

0.0 

fiplH        .Slnnp  rnpffiripnt  for  lop'-linpar  pniiatinn 

1.0 

tit r UrUL Y 

Enter  species-specific  parameters  for  polynomial  form 

MLr  UrULY 

and  defect  correction  equation  for  board  foot  volume 

(^o) 

estimates  (BFFDPOLY)  or  merchantable  cubic  foot  volume 
estimates  (MCr  JJr^ULi  i ). 
field  1:  Numeric  code  for  the  species  for  which  the 

equation  is  to  be  changed.  The  default  equation 

cnrinlipc  n  miiltinlipr  r^f  1  0  for  nil  Qnpi^ipc: 

oUpUilCO    a   1 11  Ul  Ll  LlilC;!    Ul    1  .\J    lUL    all  O^CVICO. 

field  2:  Intercept  term  for  polynomial  equation. 

1.0 

field  3:  Coefficient  for  linear  term  in  polynomial 

equation. 

0.0 

field  4:  Coefficient  for  quadratic  term  in  polynomial 

equation. 

0.0 

field  5:  Coefficient  for  cubic  term  in  polynomied 

equation. 

0.0 

field  6:  Coefficient  for  quartic  term  in  polynomial 

equation. 

0.0 

BFVOLUME 

Redefine  the  merchantability  limits  for  the  merchantable 

VOLUME 

cubic  foot  volume  equation. 

(22.  S-18) 

field  1:  Cycle  in  which  limits  defined  below 

Implement  at  start 

will  be  implemented. 

of  projection 

field  2:  Numeric  code  for  the  species  for  which 

limits  are  to  be  changed. 

Change  for  all  species 

field  3:  Minimum  d.b.h. 

6  inches  for  lodgepole  pine 

7  inches  for  all  other  species 

field  4:  Minimum  top  diameter. 

4.5  inches 

field  5:  Stump  height. 

1  foot 

MCFDLN 

Parameters  same  as  for  BFFDLN. 

MCFDPOLY 

Parameters  same  as  for  BFFDPOLY. 

VOLUME 

Parameters  same  as  for  BFVOLUME. 

31 


Specifying  Management  Activities 

Keyword 

(page  Default  parameter 

reference)  Keyword  use  and  associated  parameters  or  conditions 


CUTEFF 
(21,  S-15) 


MINHARV 
(22) 


SPECPREF 
(26) 


TCONDMLT 
(27) 


THINABA 
THINATA 
(27,  S-16) 


THINAUTO 
(28) 


Change  the  assumed  effectiveness  of  thinning  for  £ill 
thinning  activities, 
field  1:  New  value  for  global  cutting  efficiency 
parameter. 

Specify  minimum  acceptable  harvest  standards  for  board 
foot  volume,  merchantable  cubic  foot  volume,  or  basal 
area  per  acre  by  cycle, 
field  1:  The  cycle  in  which  minimum  harvest  standards 

will  be  applied, 
field  2:  The  minimum  acceptable  harvest  volume  in 

merchantable  cubic  feet  per  acre, 
field  3:  The  minimum  acceptable  heirvest  volume  in 

board  feet  per  acre, 
field  4:  The  minimum  acceptable  harvest  in  square  feet 
of  basal  area  per  acre. 
Change  the  species  component  of  the  removal  priority 
formula. 

field  1:  Date  at  which  change  is  to  be  implemented. 

field  2:  Numeric  code  for  species  whose  removal  priority 
is  to  be  changed. 

field  3:  Species  preference  vsdue. 
Change  the  impact  of  tree  value  class  on  the  determination 
of  removal  priority. 

field  1:  Date  at  which  change  is  to  be  implemented. 

field  2:  New  tree  condition  class  multiplier. 
Schedule  thinning  from  above  to  a  basal  area  per  acre 
(THINABA)  or  a  trees  per  acre  (THINATA)  target. 

field  1:  Date  that  thinning  is  scheduled. 

field  2:  The  residual  stand  density. 

field  3:  Cutting  efficiency  parameter  specific  to 

this  thinning  request, 
field  4:  d.b.h.  of  smallest  tree  that  will  be  cut. 
field  5:  d.b.h.  of  largest  tree  that  will  be  cut. 

Schedule  automatic  stocking  control.  As  nearly  as  is 
possible,  stand  density  will  be  maintained  within  a  range 
determined  by  the  minimum  and  maximum  percentage  of 
normal  stocking  entered  in  fields  2  and  3. 
field  1:  Date  that  automatic  stocking  control  is 

scheduled  to  begin, 
field  2:  Percentage  of  normal  stocking  that  defines 

the  lower  limit  for  stand  density, 
field  3:  Percentage  of  normal  stocking  that  defines 

the  upper  limit  for  stand  density, 
field  4:  Cutting  efficiency  parameter  specific  to 
automatic  stocking  control  request. 


0.98 


Applied  in  all  cycles 
0  ft3/acre 
0  bd  ft/acre 
0  ft2/acre 


Implement  at  start 
of  projection 

Ignore  the  SPECPREF 

request 

0 


Implement  at  start 
of  projection 
100 


Schedule  at  start 
of  projection 
Ignore  the  thinning 
request 

0.98 

0  inches 
999  inches 


Begin  at  start 
of  projection 

45% 

60% 

0.98 


(con.) 
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Specifying  Management  Activities  (Con.) 

Keyword 


{page 

Default  parameter 

reference) 

Keyword  use  and  associated  parameters 

or  conditions 

THINBBA 

Schedule  thinning  from  below  to  a  basal  area  per  acre 

THINBTA 

(THINBBA)  or  trees  per  acre  (THINBTA)  target. 

(27.  S-16) 

field  1:  Date  that  thinning  is  scheduled. 

Scheduled  at  start 
of  projection 

• 

field  2:  The  residual  stand  density. 

Ignore  the  thinning 
request 

field  3:  Cutting  efficiency  parameter  specific  to 

tnib  tllllllllll^  rt;i|Uc:ot. 

field  4:  d.b.h.  of  smallest  tree  that  will  be  cut. 

0  inches 

iieiu  0.  Q.D.n.  oi  icugesL  tree  tnaL'  win  oe  cul. 

yyy  incnes 

THINDBH 

Schedule  the  removal  of  a  segment  of  the  d.b.h. 

U4) 

distribution. 

iieiQ  1.  uaLG  tnac  tmnning  is  scneuuieci. 

ocneuuiea  at  stare 
of  projection 

field  2:  Smsdlest  d.b.h.  to  be  removed. 

0  inches 

field  3:  Largest  d.b.h.  to  be  removed. 

QQQ  in/«liAO 

yyy  mcnes 

field  4:  Cutting  efficiency  parameter  specific  to 

this  request. 

0.98 

THINPRSC 

Schedule  prescription  thinning.  Harvest  trees  that  were 

(24) 

marked  for  removed  on  the  input  tree  records. 

field  1:  Date  that  prescription  thinning  is  scheduled. 

Scheduled  at  start 
of  projection 

field  2:  Cutting  efficiency  parameter  specific  to 

this  thinning  request. 

0.98 

Controlling  Program  Output 

Keyword 

(page 

Default  parameter 

reference) 

Keyword  use  and  associated  parameters 

or  conditions 

COMMENT 

Enter  a  comment  that  will  be  reproduced  in  the  Input 

(48) 

Summary  Table. 
Supplemental  records:  Enter  your  comment  using  all  80 
columns  on  as  many  records  as  desired.  Signify  the 
end  of  your  comment  by  supplying  a  record  with  the 
word  "END"  entered  in  the  first  3  columns.  The  4th 

column  must  be  blank. 

None 

ECHOSUM 

Request  that  summary  output  be  copied  to  a  retrievable 

(48) 

data  storage  file. 

field  1:  Dataset  reference  number  for  output  file. 

4 

TREELIST 

Print  a  list  of  all  sample  tree  records. 

(47.  S-14) 

field  1:  Cycle  in  which  tree  list  is  to  be  printed. 

Print  tree  list  in 
all  cycles 

field  2:  Dataset  reference  number. 

3 

field  3:  If  any  number  is  entered  in  field  3,  tree 

list  will  be  printed  without  headings. 

Print  headings 
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Linkage  to  Prognosis  Model  Extensions 


Keyword 

(page 
reference) 


Keyword  use  and  associated  parameters 


Default  parameter 
or  conditions 


CHEAPO 
(86) 


COVER 
(86,  S-16) 


DFTM 
(85) 
END 
(85) 

ESTAB 

(86) 

MPB 

(85) 

WSBW 
(85) 


Generate  output  file  required  for  subsequent  execution  of 
the  CHEAPO  economic  analysis  program, 
field  1:  Dataset  reference  number  for  CHEAPO  output 
file. 

Indicates  the  start  of  special  keyword  input  file  for  the 
CANOPY  and/or  SHRUBS  extensions. 

field  1:  Cycle  to  begin  COVER  calculations. 

field  2:  Dataset  reference  number  for  COVER  output. 
Indicates  start  of  special  keyword  input  file  for  the 
Douglas-fir  tussock  moth  extension. 
Indicates  end  of  special  keyword  input  fUe  for  any 
extension. 

Indicates  start  of  special  keyword  input  file  for  the 
regeneration  establishment  extension. 
Indicates  start  of  special  keyword  input  file  for  the 
mountain  pine  beetle  extension. 

Indicates  start  of  special  keyword  input  file  for  the  western 
spruce  budworm  extension. 


11 


Beginning  of  projection 

18 


Growth  Prediction  Modifiers  and  Special  Input/Output  Options 


Keyword 

(page 
reference) 


Keyword  use  and  associated  parameters 


Default  parameter 
or  conditions 


ADDFILE 
(95) 

BAIMULT 

HTGMULT 

REGDMULT 

REGHMULT 

(94) 


BAMAX 
(95) 


COMPRESS 
(S-19) 


Specify  a  dataset  reference  number  for  a  supplemental 
keyword  record  file. 

field  1:  Dataset  reference  number. 
Enter  multiplier  to  change  prediction  of  tree  basal  area 
increment  (BAIMULT),  large  tree  height  increment 
(HTGMULT),  small  tree  diameter  increment  (REGDMULT), 
or  small  tree  height  increment  (REGHMULT). 

field  1:  Cycle  in  which  growth  multiplier  takes  effect. 

field  2:  Numeric  code  for  species  to  which  growth 
multiplier  is  to  be  applied. 

field  3:  Growth  multiplier. 
Modify  the  maximum  basal  area  used  to  control  mortality 
predictions. 

field  1:  Maximum  basal  area. 

Reduce  the  number  of  tree  records  used  to  represent  the 
stand. 

field  1:  Cycle  in  which  tree  list  will  be  compressed, 
field  2:  Number  of  tree  records  that  will  remain  after 
compression. 

field  3:  Percentage  of  new  records  that  will  be  determined 
by  finding  the  largest  gaps  in  the  classification 
space;  remainder  will  be  determined  by  splitting 
the  classes  with  the  most  variation. 

field  4:  Numeric  entry  prints  debug. 

Output  from  compression  algorithm. 


None 


Apply  in  all  cycles 

Apply  to  all  species 
1.0 


See  table  17,  Wykoff 
and  others  (1982) 


Beginning  of  projection 


150 


50 


No  Debug  output 


(con. 
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Growth  Prediction  Modifiers  and  Special  Input/Output  Options  (Con.) 

Keyword 

(page  Default  parameter 

reference)  Keyword  use  and  associated  parameters  or  conditions 


DATELIST  Instruct  program  to  print  date  of  last  revision  for 
(96)  Prognosis  Model  subprograms  and  common  areas. 

DEBUG  Request  printout  of  the  results  of  most  program  calculations 

(96)  in  any  or  all  cycles. 

field  1:  Cycle  in  which  debug  output  is  to  be  printed. 
DGSTDEV         Change  the  limits  of  the  Normal  distribution  from  which 
(93)  random  errors  are  drawn  for  increment  predictions. 

field  1:  Number  of  standard  deviations  that  define 
the  bounds  of  distribution. 
HTGMULT         Parameters  same  as  for  BAIMULT. 
HTGSTOP  Set  height  increment  to  zero  or  kill  top  for  a  random 

TOPKILL  sample  of  trees. 

(S- 20)  field  1:  Cycle  in  which  top-kill  or  growth  loss  will  be 

simulated, 
field  2:  Species  that  will  be  affected, 
field  3:  Lower  limit  of  range  of  heights  that  will  be 

affected. 

field  4:  Upper  limit  of  range  of  heights  that  will  be 
affected. 

field  5:  Probability  that  a  given  tree  will  receive 

top-kill  or  growth  loss, 
field  6:  Mean  of  the  distribution  of  the  proportional 

total  tree  height  lost  to  top-kill  (TOPKILL 

record  only). 

field  7:  Standard  deviation  of  the  distribution  of  the 
proportion  of  total  tree  height  lost  to  top-kill 
(TOPKILL  record  only). 
MORTMULT      Enter  multipher  to  change  mortality  rate  predictions. 
(94,  S-17)  field  1:  Cycle  in  which  multiplier  takes  effect. 

field  2:  Numeric  code  for  species  to  which  multiplier 

is  to  be  applied, 
field  3:  Mort£ility  multiplier. 

field  4:  d.b.h.  of  the  sm£illest  tree  to  which  multiplier 

will  be  applied, 
field  5:  d.b.h.  of  the  largest  tree  to  which  multiplier 

will  be  applied. 

NOCALIB  Suppress  calculation  of  scale  factors  for  large  tree 

(90)  diameter  increment  model  and  small  tree  height 

increment  model. 

NOTREES  Begin  a  projection  with  no  tree  records.  Permits  "bare 

(S-00)  ground"  simulations  with  records  generated  by  the 

Regeneration  EstabUshment  model. 
NOTRIPLE        Suppress  tree  record  tripling  feature. 
(93) 

NUMTRIP         Change  the  number  of  times  tree  records  will  be  tripled. 
(93)  field  1:  Number  of  triples. 

RANNSEED       Reseed  the  random  number  generator. 
(94,  S- 14)  field  1:  Replacement  seed. 

READCORD       Enter  multipliers  for  the  diameter  increment  model 
READCORH       (READCORD),  the  height  increment  model 
READCORR       (READCORH),  or  the  small  tree  height  increment  model 
(90)  (READCORR)  that  are  incorporated  prior  to  model 

caJibration. 


None 


Print  in  all  cycles 


2.0 


First  cycle 
All  species 

0.0 

0.0 

0.0 

0.0 

0.0 

Apply  in  all  cycles 

Apply  to  all  species 
1.0 

0.0 

999.0 

Calculate  scale  factors 

Two  tree  records  must  be 
input  or  projection  will 
not  run 
Tree  records 
tripled  twice 

2.0 

55329 


(con. 
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Growth  Prediction  Modifiers  and  Special  Input/Output  Options  (Con.) 


Keyword 

(page 
reference) 


Keyword  use  and  associated  parameters 


Default  parameter 
or  conditions 


REGDMULT 

REGHMULT 

REUSCORD 

REUSCORH 

REUSCORR 

(91) 

REWIND 
(95) 


STATS 
(S-19) 

TOPKILL 
(S-20) 


Supplemental  record  1:  Multipliers  for  white  pine, 
larch,  Douglas-fir,  grand  fir,  western  hemlock, 
western  redcedar,  lodgepole  pine,  and  Engelmann 
spruce. 

Supplemental  record  2:  Multipliers  for  subalpine  fir, 

ponderosa  pine,  and  mountain  hemlock. 
Pargmieters  same  as  for  BIAMULT. 
Parameters  same  as  for  BIAMULT. 

Use  multipliers  that  were  entered  with  a  READCORD,  a 
READCORH,  or  a  READCORR  in  a  previous  projection 
in  the  same  runstream. 

Causes  the  computer  to  move  the  read  position  pointer  to 
the  beginning  of  the  dataset  referenced  by  the  unit 
number  entered  in  field  1.  This  record  is  useful  when 
multiple  projections  are  made  with  the  same  tree  record 
file  in  a  single  runstream. 

field  1:  Dataset  reference  number. 
Prepare  and  print  table  that  describes  distribution  of 
input  data. 

field  1:  Significance  level  for  computing  confidence  limits. 
Parameters  same  as  for  HTGSTOP. 


Default  value  for 
£dl  multipliers  is  1.0 


0.05 
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Wykoff,  William  R.  Supplement  to  the  user's  guide  for  the  Stand  Prognosis 
Model— version  5.0.  General  Technical  Report  INT-208.  Ogden,  UT;  U.S. 
Department  of  Agriculture,  Forest  Service,  Intermountain  Research  Station; 
1986.  36  p. 

Differences  between  Prognosis  Model  versions  4.0  and  5.0  are  described. 
Additions  to  version  5.0  include  an  event  monitor  that  schedules  activities  con- 
tingent on  stand  characteristics,  a  regeneration  establishment  model  that 
predicts  the  structure  of  the  regeneration  stand  following  treatment,  and  a 
COVER  model  that  predicts  shrub  development  and  total  canopy  cover.  Program 
performance  has  been  enhanced  by  modifications  to  several  of  the  submodels 
that  predict  tree  increments. 


KEYWORDS:  growth  and  yield,  forest  management,  planning,  growth  projection, 
stand  models,  tree  increments,  tree  mortality 


INTERMOUNTAIN  RESEARCH  STATION 


The  Intermountain  Research  Station  provides  scientific  knowl- 
edge and  technology  to  improve  management,  protection,  and  use 
of  the  forests  and  rangelands  of  the  Intermountain  West.  Research 
is  designed  to  meet  the  needs  of  National  Forest  managers, 
Federal  and  State  agencies,  industry,  academic  institutions,  public 
and  private  organizations,  and  individuals.  Results  of  research  are 
made  available  through  publications,  symposia,  workshops,  training 
sessions,  and  personal  contacts. 

The  Intermountain  Research  Station  territory  includes  Montana, 
Idaho,  Utah,  Nevada,  and  western  Wyoming.  Eighty-five  percent  of 
the  lands  in  the  Station  area,  about  231  million  acres,  are  classified 
as  forest  or  rangeland.  They  include  grasslands,  deserts,  shrub- 
lands,  alpine  areas,  and  forests.  They  provide  fiber  for  forest  in- 
dustries, minerals  and  fossil  fuels  for  energy  and  industrial  develop- 
ment, water  for  domestic  and  industrial  consumption,  forage  for 
livestock  and  wildlife,  and  recreation  opportunities  for  millions  of 
visitors. 

Several  Station  units  conduct  research  in  additional  western 
States,  or  have  missions  that  are  national  or  international  in  scope. 
Station  laboratories  are  located  in: 

Boise,  Idaho 

Bozeman,  Montana  (in  cooperation  with  Montana  State 
University) 

Logan,  Utah  (in  cooperation  with  Utah  State  University) 

Missoula,  Montana  (in  cooperation  with  the  University  of 
Montana) 

Moscow,  Idaho  (in  cooperation  with  the  University  of  Idaho) 
Ogden,  Utah 

Provo,  Utah  (in  cooperation  with  Brigham  Young  University) 
Reno,  Nevada  (in  cooperation  with  the  University  of  Nevada) 


